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We  present  the  measurement  of  squeezed  light  generation  using  an  engineered  optomechanical 
system  fabricated  from  a  silicon  microchip  and  composed  of  a  micromechanical  resonator  coupled  to  a 
nanophotonic  cavity.  Laser  light  is  used  to  measure  the  fluctuations  in  the  position  of  the  mechanical 
resonator  at  a  measurement  rate  comparable  to  the  free  dynamics  of  the  mechanical  resonator,  and 
greater  than  its  thermal  decoherence  rate.  By  approaching  the  strong  continuous  measurement 
regime  we  observe,  through  homodyne  detection,  non-trivial  modifications  of  the  reflected  light’s 
vacuum  fluctuation  spectrum.  In  spite  of  the  mechanical  resonator’s  highly  excited  thermal  state 
(10,000  phonons),  we  observe  squeezing  at  the  level  of  4.5  d=  0.5%  below  that  of  shot-noise  over 
a  few  MHz  bandwidth  around  the  mechanical  resonance  frequency  of  28  MHz.  This  squeezing  is 
interpreted  as  an  unambiguous  quantum  signature  of  radiation  pressure  shot-noise. 


Monitoring  a  mechanical  object’s  motion,  even  with  the 
gentle  touch  of  light,  fundamentally  alters  its  dynamics. 
The  experimental  manifestation  of  this  basic  principle  of 
quantum  mechanics,  its  link  to  the  quantum  nature  of 
light,  and  the  extension  of  quantum  measurement  to  the 
macroscopic  realm  have  all  received  extensive  attention 
over  the  last  half  century  [1,  2].  The  use  of  squeezed 
light,  with  quantum  fluctuations  below  that  of  the  vac¬ 
uum  field,  was  proposed  nearly  three  decades  ago  [3]  as  a 
means  for  beating  the  standard  quantum  limits  in  preci¬ 
sion  displacement  and  force  measurements.  Conversely,  it 
has  also  been  proposed  that  a  strong  continuous  measure¬ 
ment  of  a  mirror’s  position  with  light,  may  itself  give  rise 
to  squeezed  light  [4,  5].  In  this  Letter,  we  present  such  a 
continuous  position  measurement  using  an  engineered  op¬ 
tomechanical  system  fabricated  from  a  silicon  microchip 
and  composed  of  a  micromechanical  resonator  coupled  to 
a  nanophotonic  cavity.  Laser  light  is  used  to  measure  the 
fluctuations  in  the  position  of  the  mechanical  resonator 
at  a  measurement  rate  comparable  to  the  free  dynamics 
of  the  mechanical  resonator,  and  greater  than  its  thermal 
decoherence  rate.  By  approaching  the  strong  continuous 
measurement  regime  we  observe,  through  homo  dyne  de¬ 
tection,  non-trivial  modifications  of  the  reflected  light’s 
vacuum  fluctuation  spectrum.  In  spite  of  the  mechanical 
resonator’s  highly  excited  thermal  state  (104  phonons), 
we  observe  squeezing  at  the  level  of  4.5  ±0.5%  below  that 
of  shot-noise  over  a  few  MHz  bandwidth  around  the  me¬ 
chanical  resonance  frequency  of  ujrn/27r  ~  28  MHz.  This 
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squeezing  is  interpreted  as  an  unambiguous  quantum  sig¬ 
nature  of  radiation  pressure  shot-noise  [6]. 

The  generation  of  states  of  light  with  fluctuations  be¬ 
low  that  of  vacuum  has  been  of  great  theoretical  inter¬ 
est  since  the  1970s  [3,  7-9].  Early  experimental  work 
demonstrated  squeezing  of  a  few  percent  in  a  large  vari¬ 
ety  of  different  nonlinear  systems,  such  as  neutral  atoms 
in  a  cavity  [10],  optical  fibers  [11],  and  crystals  with  bulk 
optical  nonlinearities  [12,  13].  This  initial  research  was 
mainly  pursued  as  a  strategy  to  mitigate  the  effects  of 
shot-noise  given  the  possibility  of  improved  optical  com¬ 
munication  [7]  and  better  sensitivity  in  gravitational  wave 
detectors  [3,  8].  In  recent  years,  in  addition  to  being  in¬ 
stalled  in  gravitational  wave  detectors  [14],  squeezed  light 
has  enhanced  metrology  in  more  applied  settings  [15]. 

The  vacuum  fluctuations  arising  from  the  quantum  na¬ 
ture  of  light  determine  our  ability  to  optically  resolve 
mechanical  motion  and  set  limits  on  the  perturbation 
caused  by  the  act  of  measurement  [16].  A  well-suited 
system  to  experimentally  study  quantum  measurement 
is  that  of  cavity-optomechanics,  where  an  optical  cav¬ 
ity’s  resonance  frequency  can  be  designed  to  be  sensi¬ 
tive  to  the  position  of  a  mechanical  system.  By  moni¬ 
toring  the  phase  and  intensity  of  the  reflected  light  from 
such  a  cavity,  a  continuous  measurement  of  mechanical 
displacement  can  be  made.  Systems  operating  on  this 
simple  premise  have  been  realized  in  a  variety  of  experi¬ 
mental  settings,  such  as  in  large-scale  laser  gravitational 
wave  interferometers  [17],  microwave  circuits  with  elec¬ 
tromechanical  elements  [18,  19],  mechanical  elements  in¬ 
tegrated  with  or  comprising  Fabry-Perot  cavities  [20-22] , 
and  on-chip  nanophotonic  cavities  sensitive  to  mechani¬ 
cal  deformations  [23,  24].  Conceptually,  the  same  basic 
optomechanical  interaction  appears  when  the  collective 
motion  of  an  ultracold  gas  of  atoms  shifts  the  resonance 
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frequency  of  an  optical  cavity,  and  recent  gas-phase  ex¬ 
periments  have  shown  a  variety  of  optomechanical  effects, 
including  squeezed  light  generation  [25]. 

The  simplest  cavity-optomechanical  system  consists  of 
two  modes,  an  optical  and  a  mechanical  resonance,  and 
is  parameterized  by  their  respective  resonance  frequencies 
( uj0 ,  cjm)  and  quality  factors  (Q0,  Om),  and  the  coupling 
rate  of  intracavity  light  to  the  resonant  mechanical  mo¬ 
tion  (go).  The  Hamiltonian  describing  the  interaction  be¬ 
tween  light  and  mechanics  is  H[nt  =  hgob)ax /xzpf,  with 
the  amplitude  of  the  mechanical  resonator  motion  being 
given  by  x  =  xzpf  •  (tf  +  6),  where  xzpf  =  y/h/2mefiUJrn  is 
the  zero-point  fluctuation  and  me &  the  effective  motional 
mass  of  the  resonator.  Here  a  (a*-)  and  b  (£d)  are  the 
annihilation  (creation)  operators  of  optical  and  mechan¬ 
ical  excitations,  respectively.  The  optical  cavity  decay 
rate,  k  =  uj0/Q0,  is  the  loss  rate  of  photons  from  the 
cavity  and  the  rate  at  which  optical  vacuum  fluctuations, 
or  shot-noise,  is  coupled  into  the  optical  resonance  [26]. 
Similarly,  the  mechanical  damping  rate  71  =  oom/Qm  is 
the  rate  at  which  thermal  bath  fluctuations  couple  into 
the  mechanical  system.  In  all  experimental  realizations  of 
optomechanics  to  date,  including  that  presented  here,  the 
optomechanical  coupling  rate  go  has  been  much  smaller 
than  the  cavity  decay  rate  k.  As  such,  without  a  strong 
coherent  drive,  the  interaction  of  the  vacuum  fluctuations 
with  the  mechanics  is  negligible. 

Under  the  effect  of  a  coherent  laser  drive,  the  cavity  is 
populated  with  a  mean  intracavity  photon  number  (nc), 
and  one  considers  the  optical  fluctuations  about  the  clas¬ 
sical  steady-state,  a  y/ (nc)  +  a.  This  modifies  the 
optomechanical  interaction  resulting  in  a  linear  coupling 
between  the  fluctuations  of  the  intracavity  optical  field 
(Xq  =  a  +  at)  and  the  position  fluctuations  of  the  me¬ 
chanical  system  x:  H[nt  =  HGX o%/xzpf.  The  paramet¬ 
ric  linear  coupling  occurs  at  an  effective  interaction  rate 
of  G  =  yj (nc)go-  Through  this  interaction,  the  inten¬ 
sity  fluctuations  of  the  vacuum  field  XQ=0(t)  entering  the 
cavity  impart  a  force  on  the  mechanical  system, 

FBA(t)  =  h'Vr™*X™(t),  (1) 

^zpf 

which  is  the  radiation  pressure  shot-noise  (RPSN)  force. 
The  mechanical  motion  in  turn  is  recorded  in  the  phase 
of  the  light  leaving  the  cavity, 

Mout\t)  =  -x£n\t)  -  2 •  sin (6),  (2) 

^zpf 

where  6  is  the  quadrature  angle  of  the  detected  optical 
field,  with  9  =  0  (tt/2)  referring  to  the  intensity  (phase) 
quadrature.  The  latter  can  be  interpreted  as  the  infor¬ 
mational  aspect  of  measurement,  with  the  optical  cavity 
playing  the  role  of  the  position  detector  (measuring  ob¬ 
servable  x  at  a  rate  Tmeas  =  4 G2/k),  while  the  former  is 
the  quantum  measurement  back-action  imposing  a  fluc¬ 
tuating  noise  force  onto  the  mechanical  system  [2] .  In  ad¬ 
dition  to  the  back-action  noise,  thermal  fluctuations  from 


the  bath  also  drive  the  mechanical  motion,  with  their 
magnitudes  becoming  comparable  as  Tmeas  approaches 
the  thermalization  rate  Tthermai  (c<;)  =  jin(u))-  The  bath 
occupation  number  n( uS)  is  equivalent  to  its  high  temper¬ 
ature  limit  in  these  measurements. 

The  auto-correlation  of  the  optical  output  field  has  con¬ 
tributions  due  to  the  auto-correlation  of  the  optical  in¬ 
put  field  and  the  mechanical  fluctuations,  as  well  as  their 
cross-correlations.  Due  to  the  back-action  force  on  the 
mechanical  system,  the  motion  has  non-zero  correlations 
with  the  RPSN  ((x(t)FBA(t'))),  which  can  be  negative 
for  some  values  of  0.  Under  certain  conditions  this  can 
lead  to  the  squeezing  of  the  optical  output  field.  The 
spectral  density  of  a  balanced  homodyne  measurement  of 
the  output  field,  normalized  to  shot-noise,  is  given  in  the 
quasi-static  limit  for  frequencies  uj  <C  c^m  (see  the  Ap¬ 
pendix)  : 

S?jV)  =  1  +  4(rmeas/wm)  sin(20) 

+4rmeasnM  (1  +  cos(26)))j  (3) 

Grri 

Measurement  of  the  spectrum  below  shot-noise  (squeez¬ 
ing)  is  therefore  attributable  to  cross-correlations  of  mea¬ 
surement  back-action  noise  due  the  optical  vacuum  fluc¬ 
tuations  and  position  fluctuations  of  the  mechanical  res¬ 
onator. 

The  primary  technical  hurdle  to  observing  such  squeez¬ 
ing,  as  in  many  quantum  measurements,  is  the  strong  cou¬ 
pling  of  a  preferred  detection  channel  (the  optical  probe) 
simultaneous  with  the  minimization  of  unwanted  environ¬ 
mental  perturbations  of  the  mechanical  system.  As  indi¬ 
cated  in  Eq.  (11),  fluctuations  from  the  thermal  bath  limit 
squeezing  of  the  optical  probe  field  to  a  regime  in  which 
n(uj)  <  Qm.  This  requirement  is  equivalent  to  having  a 
Q-frequency  product,  Qm^m  >  T^/h.  Equation  (11) 

also  indicates  that  squeezing  becomes  appreciable  only  as 
the  measurement  rate  of  the  mechanics  by  the  light  field 
approaches  the  mechanical  frequency.  In  some  respects, 
both  of  these  challenges  have  been  overcome  in  recent  op¬ 
tomechanics  work  [19,  24,  27],  whereby  thermal  bath  cou¬ 
pling  has  been  made  smaller  than  the  optically-induced 
cooling  of  the  mechanical  resonator.  However,  significant 
squeezing  over  an  appreciable  spectral  bandwidth  requires 
not  only  a  large  cooperativity  between  light  and  mechan¬ 
ics,  as  represented  by  C  =  Tmeas/7 i  and  realized  in  cool¬ 
ing  experiments,  but  the  more  stringent  requirement  that 
the  effective  measurement  back-action  force  be  compa¬ 
rable  to  all  forces  acting  on  the  mechanics,  including  the 
elastic  restoring  force  of  the  mechanical  structure.  A  more 
comprehensive  model  (see  Appendix),  including  the  reso¬ 
nant  response  of  the  mechanical  system  to  RPSN,  shows 
that  the  strict  requirements  of  the  quasi-static  model  are 
somewhat  relaxed,  and  that  the  squeezing  scales  approx¬ 
imately  as  rmeas/£,  where  5  is  the  effective  bandwidth  of 
squeezing  around  the  mechanical  resonance  frequency. 

In  order  to  meet  the  requirements  of  strong  measure¬ 
ment  and  efficient  detection,  we  designed  a  zipper-style 
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FIG.  1:  Optomechanical  device,  a,  Scanning  electron  mi¬ 
croscope  image  of  a  waveguide-coupled  zipper  optomechanical 
cavity.  The  waveguide  width  is  adiabatically  tapered  along  its 
length  and  terminated  with  a  photonic  crystal  mirror  next  to 
the  cavity.  The  tapering  of  the  waveguide  allows  for  efficient 
input/output  coupling  while  the  photonic  crystal  termination 
makes  the  coupling  to  the  cavity  single-sided.  Two  zipper  cav¬ 
ities  are  coupled  above  and  below  the  waveguide,  each  with  a 
slightly  different  optical  resonance  frequency  allowing  them  to 
be  separately  addressed,  b,  (left)  Close-up  of  the  coupling  re¬ 
gion  between  one  of  the  cavities  and  the  waveguide,  (right) 
Finite  element  method  (FEM)  simulation  of  the  cavity  field 
leaking  into  the  waveguide  (log  scale).  Note  that  the  field 
does  not  leak  into  the  mirror  region  of  the  waveguide,  c,  (top) 
FEM  simulation  showing  the  in-plane  electrical  field  of  the 
fundamental  optical  cavity  mode,  (bottom)  FEM  simulation 
of  the  displacement  of  the  fundamental  in-plane  differential 
mode  of  the  structure  with  frequency  um/27i  =  28  MHz.  The 
mechanical  motion,  modifying  the  gap  between  the  beams, 
shifts  the  optical  cavity  frequency  leading  to  optomechanical 
coupling. 


optomechanical  cavity  [23]  with  a  novel  integrated  waveg¬ 
uide  coupler  fabricated  from  the  220  nm  thick  silicon  de¬ 
vice  layer  of  a  silicon-on-insulator  microchip  (see  Fig.  la). 
The  in-plane  differential  motion  of  the  two  beams  at  a 
fundamental  frequency  of  c<;m/27r  =  28  MHz  strongly 
modulates  the  co-localized  fundamental  optical  resonance 
of  the  cavity  with  a  theoretical  vacuum  coupling  rate  of 
go/2ir  =  1  MHz.  As  shown  in  Fig.  lb,  we  use  a  silicon 
waveguide  with  a  high  reflectivity  photonic  crystal  end- 
mirror  to  efficiently  excite  and  collect  light  from  the  op¬ 
tical  cavity.  Coupling  light  from  the  silicon  waveguide  to 
a  single-mode  optical  fiber  is  performed  using  an  optical 
fiber  taper  and  a  combination  of  adiabatic  mode  coupling 
and  transformation  (see  Fig.  lb). 

The  experimental  setup  used  to  characterize  the  zipper 
cavity  system  and  measure  the  optomechanical  squeezing 
of  light  is  shown  in  Fig.  2a.  The  silicon  sample  is  placed 
in  a  continuous  flow  4He  cryostat  with  a  cold  finger  tem¬ 
perature  of  10  K.  A  signal  laser  beam  is  used  to  probe 
the  optomechanical  system  and  measure  the  mechanical 


motion  of  the  zipper  cavity.  A  wavelength  scan  of  the  re¬ 
flected  signal  from  the  cavity  is  plotted  in  Fig.  2c,  showing 
an  optical  resonance  with  a  linewidth  k/2tt  =  3.42  GHz 
at  a  wavelength  of  Ac  =  1540  nm.  Inefficiencies  in  the 
collection  and  detection  of  light  correspond  to  additional 
uncorrelated  shot-noise  in  the  signal  and  can  reduce  the 
squeezing  to  undetectable  levels.  For  the  device  studied 
here,  the  cavity  coupling  efficiency,  corresponding  to  the 
percentage  of  photons  sent  into  the  cavity  which  are  re¬ 
flected,  is  determined  to  be  rjk  =  0.54.  The  fiber-to-chip 
coupling  efficiency  is  measured  at  rjcp  =  0.90.  A  homo¬ 
dyne  detection  scheme  [28]  allows  for  high  efficiency  de¬ 
tection  of  arbitrary  quadratures  of  the  optical  signal  field. 
Characterization  and  optimization  of  the  efficiency  of  the 
entire  optical  signal  path  and  homodyne  detection  sys¬ 
tem  (see  Appendix  for  details)  results  in  an  overall  setup 
efficiency  of  r]setup  =  0.48,  corresponding  to  a  total  signal 
detection  efficiency  of  ^ot  =  VsetupV^  =  0.26. 

Figure  2c  shows  the  noise  spectrum  of  the  thermal  mo¬ 
tion  of  the  mechanical  resonator  obtained  by  positioning 
the  laser  frequency  near  the  cavity  resonance  and  tuning 
the  relative  local  oscillator  (LO)  phase  of  the  homodyne 
detector,  #iock,  to  measure  the  quadrature  of  the  reflected 
signal  in  which  mechanical  motion  is  imprinted  (roughly 
the  phase  quadrature  for  near-resonance  probing).  The 
mechanical  spectrum  is  seen  to  contain  the  in-plane  dif¬ 
ferential  mode  of  interest  at  cjm/27r  =  28  MHz,  as  well 
as  several  other  more  weakly  coupled  mechanical  reso¬ 
nances  of  the  nanobeams  and  coupling  waveguide  (the 
in-plane  differential  mode  peak  appears  reduced  relative 
to  the  other  modes  in  this  plot  due  to  the  limited  resolu¬ 
tion  bandwidth  of  the  measurement).  A  high-resolution, 
narrowband  spectrum  of  the  in-plane  differential  mode  is 
displayed  as  an  inset  to  Fig.  2c,  and  shows  a  linewidth  of 
7^/27t  =  172  Hz,  corresponding  to  a  mechanical  Q-factor 
of  Qm  =  1.66  x  105.  The  vacuum  coupling  rate  of  the 
in-plane  differential  mode,  measured  from  the  detuning 
dependence  of  the  optical  spring  shift  and  damping  (see 
Appendix),  is  determined  to  be  go/2ir  =  750  kHz,  in  good 
correspondence  with  theory.  From  the  calibration  of  the 
noise  power  under  the  Lorentzian  in  Fig.  2c,  the  in-plane 
differential  mode  is  found  to  thermalize  (at  low  optical 
probe  power)  to  a  temperature  of  ~  16  K,  correspond¬ 
ing  to  a  phonon  occupancy  of  (n)  ~  1.2  x  104.  This  yields 
a  ratio,  Qm^m/fe^b  ~  13,  well  within  the  regime  where 
coherent  motion  and  squeezing  are  possible. 

In  order  to  systematically  and  accurately  study  the 
noise  properties  of  the  reflected  optical  signal  from  the 
cavity  we  make  a  series  of  measurements  to  character¬ 
ize  our  laser  and  detection  setup.  Figure  2d  shows  the 
measured  noise  power  spectral  density  (PSD)  of  the  bal¬ 
anced  homodyne  detector  (dark  current  subtracted)  for 
uj  ~  as  a  function  of  LO  power  (signal  blocked), 
indicating  a  linear  dependence  on  power  and  negligible 
added  noise  above  shot-noise.  In  the  measured  squeezing 
data  to  follow,  a  LO  power  of  3  mW  is  used.  Calibra¬ 
tion  of  the  laser  intensity  and  frequency  noise  over  the 
frequency  range  of  interest  (uj/2tt  =  1-40  MHz)  is  mea- 
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FIG.  2:  Experimental  setup  and  device  characterization,  a  A  tunable  external-cavity  diode  laser,  actively  locked  to 
a  wavemeter  (A-meter),  is  used  to  generate  a  strong  local  oscillator  (LO)  and  the  measurement  signal.  A  portion  of  each 
field  is  split  off,  detected,  and  used  to  individually  stabilize  the  intensity  in  each  arm  using  intensity  modulators  (IM).  Fiber 
polarization  controllers  (FPC)  adjust  the  polarization  of  the  LO  and  signal.  A  variable  optical  attenuator  (VOA)  is  used  to 
set  the  signal  power  and  an  acousto-optic  modulator  (AOM)  is  used  to  generate  a  tone  of  known  amplitude  for  calibration  (see 
Appendix).  The  signal  beam  is  fed  into  a  fiber  taper  that  is  mounted  inside  a  4He  continuous  flow  cryostat  and  is  coupled  to 
the  zipper  cavity.  The  reflected  signal  beam  from  the  cavity  is  separated  using  a  circulator,  and  is  switched  between  one  of 
three  detection  paths:  a  power  meter  (PM)  for  power  calibration,  a  photodetector  (PD1)  for  optical  spectroscopy  of  the  cavity, 
and  an  erbium-doped  fiber  amplifier  (EDFA)  for  spectroscopy  of  the  coupled  cavity-optomechanical  system  on  a  real-time 
spectrum  analyzer  (RSA)  or  network  analyzer  (ENA).  For  the  measurement  of  squeezing,  the  reflected  cavity  signal  is  sent  to 
a  fourth  path  containing  a  variable  coupler  (VC),  where  it  is  recombined  with  the  strong  local  oscillator  and  detected  on  a 
balanced  homodyne  detector  (BHD).  The  relative  phase  between  the  LO  and  the  signal  is  set  using  a  fiber  stretcher  (FS).  b 
(top)  Reflected  signal  from  the  optical  cavity  at  low  optical  power  (nc  ~  10,  linewidth  k,/2tt  =  3.42  GHz),  (bottom)  High- 
power  (nc  =  790)  reflected  signal,  showing  the  cavity-laser  detuning  (dashed  line)  locked  to  during  squeezing  measurements,  c 
Homodyne  noise  PSD  of  the  reflected  signal  showing  the  transduced  thermal  Brownian  motion  of  the  zipper  cavity  at  Tb  =  16  K 
(green  curve;  nc  =  80).  The  red  curve  is  the  shot  noise  level  and  the  black  curve  is  the  detector’s  dark  noise.  The  inset  shows 
a  zoom-in  of  the  fundamental  in-plane  differential  mechanical  mode  of  the  zipper  cavity  (linewidth  7i/27r  =  172  Hz),  d  Mean 
value  of  the  PSD  of  the  balanced  homodyne  detector  as  a  function  of  the  LO  power  (signal  blocked).  The  filled  data  point 
indicates  the  LO  power  used  in  the  squeezing  measurements.  The  red  and  dashed  black  curves  correspond  to  a  linear  fit  to 
the  data  and  the  detector  dark  current  level,  respectively,  e  Noise  PSD  as  a  function  of  6\oc k  with  the  signal  detuned  far 
off-resonance  at  A/k  ~  30  (shades  green  to  red),  referenced  to  the  noise  level  with  the  signal  blocked  (blue). 


sured  by  direct  photodetection  of  the  laser,  pre-  and  post- 
transmission  through  a  fiber  Mach-Zehnder  interferome¬ 
ter  with  a  known  frequency  response.  The  laser  intensity 
noise  is  measured  to  be  shot-noise  dominated  over  this  fre¬ 
quency  range,  while  the  laser  frequency  noise  is  measured 
to  be  roughly  flat  at  a  level  of  SULj  ^  5  x  103  rad2 -Hz. 
Laser  phase  noise  on  the  signal  beam  can  be  converted 
to  intensity  noise  by  reflection  from  the  dispersive  cavity 
or  due  to  frequency  dependent  components  in  the  optical 
train,  and  can  add  to  the  detected  noise  floor.  Due  to 
the  broad  linewidth  of  the  zipper  cavity  resonance,  the 
effects  of  the  laser  phase  noise  are  found  to  be  negligi¬ 
ble.  The  full  suite  of  noise  and  calibration  measurements 


performed  are  described  in  detail  in  the  Appendix. 

Measurement  of  the  noise  in  the  reflected  optical  signal 
from  the  cavity  as  a  function  of  quadrature  angle,  fre¬ 
quency,  and  signal  power  is  presented  in  Figs.  3  and  4. 
These  measurements  are  performed  for  laser  light  on  res¬ 
onance  with  the  optical  cavity  and  for  input  signal  powers 
varying  from  252  nW  to  3.99  fi W  in  steps  of  2  dB,  with  the 
maximum  signal  power  corresponding  to  an  average  intra¬ 
cavity  photon  number  of  (nc)  =  3, 153.  Positioning  of  the 
laser  at  the  appropriate  cavity  detuning  for  each  signal 
power  is  performed  by  scanning  the  wavelength  across 
the  cavity  resonance  while  monitoring  the  reflection,  and 
then  stepping  the  laser  frequency  towards  the  cavity  from 
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FIG.  3:  Optomechanical  squeezing  of  light,  a,  Theoretical  model.  Density  plot  of  the  predicted  reflected  signal  noise 
PSD,  as  measured  on  a  balanced  homodyne  detector  and  normalized  to  shot-noise,  for  a  simplified  model  of  the  optomechanical 
system  (see  Appendix).  Areas  below  shot-noise  are  shown  in  blue  shades  on  a  linear  scale.  Areas  with  noise  above  shot-noise 
are  shown  in  orange  shades  on  a  log-scale.  The  solid  white  line  is  a  contour  delineating  noise  above  and  below  shot-noise,  b, 
Experimental  data.  Density  plot  of  the  measured  reflected  signal  noise  PSD  for  nc  —  790  normalized  to  the  measured  shot-noise 
level,  c,  Slice  of  the  measured  density  plot  in  b  taken  at  ^lock/^r  =  0.23.  d,  Slice  of  the  measured  density  plot  in  b  taken  at 
#iock/7r  =  0.16.  In  c  and  d,  the  black  curve  corresponds  to  the  measured  data  slice  extracted  from  b.  The  dark  blue  traces  are 
several  measurements  of  the  shot-noise  level  (average  shown  in  light  blue).  Also  indicated  is  a  model  of  the  squeezing  in  the 
absence  of  thermal  noise  (red  dashed  curve),  the  same  model  with  thermal  noise  included  (solid  red  curve),  and  a  full  noise 
model  including  additional  phenomenological  noise  sources  (solid  green  curve). 


the  red  side  until  the  reflection  matches  the  level  that  cor¬ 
responds  to  a  detuning  of  Aiock/^  ~  0.04  (see  Fig.  2b). 
The  laser  is  locked  to  this  frequency  using  a  wavemeter 
with  a  frequency  resolution  of  ±0. 0015ft.  Drift  of  the  op¬ 
tical  cavity  resonance  over  the  timescale  of  a  single  noise 
spectrum  measurement  (minutes)  is  found  to  be  negligi¬ 
ble.  An  estimate  of  the  variance  of  Aiock  is  determined 
from  the  dependence  of  the  transduction  of  the  mechani¬ 
cal  motion  on  the  quadrature  phase,  indicating  that  from 
one  lock  to  another  Aiock/ft  =  0.044  db  0.006. 

In  Fig.  3  we  plot  the  theoretically  predicted  and  mea¬ 
sured  noise  PSD  versus  quadrature  angle  for  a  signal 
power  corresponding  to  (nc)  =  790  photons.  Each 
quadrature  spectrum  is  the  average  of  150  traces  taken 
over  20  seconds,  and  after  every  other  spectrum,  the  sig¬ 
nal  arm  is  blocked  and  the  shot-noise  PSD  is  measured. 
The  shot-noise  level,  which  represents  the  noise  of  the 
electromagnetic  vacuum  on  the  signal  arm,  is  used  to  nor¬ 
malize  the  spectra.  We  find  at  certain  quadrature  angles, 
and  for  frequencies  a  few  MHz  around  the  mechanical  res¬ 
onance  frequency,  that  the  light  reflected  from  the  zipper 
cavity  shows  a  noise  PSD  below  that  of  vacuum.  The  den¬ 
sity  plot  of  the  theoretically  predicted  noise  PSD  (Fig.  3a) 
shows  the  expected  wideband  squeezing  due  to  the  strong 
optomechanical  coupling  in  these  devices,  as  well  as  a 
change  in  the  phase  angle  where  squeezing  is  observed  at 
below  and  above  the  mechanical  frequency.  This  change  is 
due  to  the  change  in  sign  of  the  mechanical  susceptibility 


and  the  corresponding  change  in  phase  of  the  mechanical 
response  to  RPSN.  The  measured  noise  PSD  density  plot 
(Fig.  3b)  shows  the  presence  of  several  other  mechanical 
noise  peaks  and  a  reduced  squeezing  bandwidth,  yet  the 
overall  phase-  and  frequency-dependent  characteristics  of 
the  squeezing  around  the  strongly-coupled  in-plane  me¬ 
chanical  mode  are  clearly  present.  In  particular,  Figs.  3c 
and  d  show  two  slices  of  the  noise  PSD  density  plot  which 
show  the  region  of  squeezing  change  from  being  below  to 
above  the  mechanical  resonance  frequency. 

In  Fig.  4a  we  show  the  measured  noise  PSD  (grey  cir¬ 
cles)  as  a  function  of  quadrature  angle  for  a  frequency 
slice  at  uj/2tt  =  27.9  MHz  of  the  data  shown  in  Fig.  3b. 
The  measured  squeezing  (anti-squeezing)  is  seen  to  be 
smaller  (larger)  than  expected  from  a  model  of  the  op¬ 
tomechanical  cavity  without  thermal  noise.  We  also  plot 
in  Fig.  4b  the  maximum  measured  and  modeled  squeez¬ 
ing  as  a  function  of  signal  power.  The  simple  theory  pre¬ 
dicts  a  squeezing  level  (blue  curve)  which  monotonically 
increases  with  signal  power,  whereas  the  measured  max¬ 
imum  squeezing  saturates  at  a  level  of  4.5  ±  0.5%  below 
the  shot-noise  at  an  intracavity  power  (nc)  =  1,984  pho¬ 
tons.  In  order  to  understand  the  noise  processes  that  limit 
the  bandwidth  and  magnitude  of  the  measured  squeez¬ 
ing,  we  plot  in  Fig.  4c  the  noise  PSD  (shot-noise  sub¬ 
tracted)  for  phase  quadratures  that  maximize  (left  plot) 
and  minimize  (right  plot)  the  transduction  of  the  me¬ 
chanical  mode  peak  at  cjm/27r  =  28  MHz.  Along  with 
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FIG.  4:  Spectral  and  power  dependence  of  noise,  a, 

Measured  (filled  circles)  balanced  homodyne  noise  power  of 
the  reflected  signal  at  uj/2t:  =  27.9  MHz  versus  quadrature 
angle  (Aiock/ft  =  0.044  and  nc  =  790).  The  red  curve  cor¬ 
responds  to  the  full  noise  model.  The  solid  blue  curve  is  for 
a  model  including  the  response  of  the  mechanical  mode  in 
the  absence  of  thermal  noise,  i.e.,  driven  by  RPSN  only  (the 
dashed  blue  curve  shows  the  thermal  noise  component),  b, 
Measured  (filled  circles)  minimum  noise  PSD  normalized  to 
shot-noise  versus  nc.  The  left  plot  is  the  maximum  squeezing 
for  uj  <  Um  and  the  right  for  uj  >  cjm.  Also  shown  is  the 
single-mode  noise  model  (blue  curve)  and  the  full  noise  model 
(red  curve),  c,  Balanced  homodyne  noise  PSD  of  the  reflected 
cavity  signal  for  Aiock/ft  =  0.052  and  nc  —  3, 153.  Left  (right) 
plot  shows  phase  quadrature  corresponding  to  maximum  (min¬ 
imum)  transduction  of  mechanical  motion.  The  black  curve  is 
the  measured  data  with  the  shot-noise  level  subtracted.  Also 
shown  are  modeled  laser  phase  noise  (green  curve) ,  the  single¬ 
mode  noise  model  (blue  curve),  and  the  full  noise  model  (red 
curve) . 


the  measured  data  (black  curve),  we  also  plot  the  esti¬ 
mated  noise  due  to  phase  noise  of  the  signal  laser  (green 
curve)  and  that  for  a  single  mechanical  mode  (blue  curve) 
assuming  a  thermal  bath  temperature  of  T&  =  16  K  and 
a  frequency-independent  damping  rate.  The  single-mode 
noise  model  greatly  underestimates  the  background  noise 
level,  especially  in  the  quadrature  minimizing  transduced 
motion  in  which  we  measure  a  a;-1/1 2  (as  opposed  to  uj~x) 


frequency  dependence  to  the  low  frequency  noise.  As  de¬ 
scribed  in  the  Appendix,  the  additional  measured  noise  is 
thought  to  arise  from  a  combination  of  the  thermal  noise 
tails  of  higher  frequency  mechanical  modes  and  fluctu¬ 
ations  in  the  optical  cavity  damping  rate,  along  with  a 
small  amount  of  heating  due  to  optical  absorption  in  the 
silicon  cavity.  The  red  curves  in  each  of  the  plots  in  Fig.  4 
show  the  full  noise  model  incorporating  these  added  phe¬ 
nomenological  terms. 

Somewhat  surprisingly,  these  measurements  show  that 
by  reflecting  light  off  a  thin-film  mechanical  resonator  un¬ 
dergoing  large  amplitude  thermal  motion,  light  that  is  in 
certain  respects  quieter  than  vacuum  can  be  obtained. 
This  is  found  to  result  from  radiation-pressure  fluctua¬ 
tions  being  strongly  imprinted  on  and  modified  by  the 
motion  of  the  mechanical  resonator,  thus  demonstrating 
a  fundamentally  quantum  aspect  of  displacement  mea¬ 
surement.  The  devices  in  this  work  utilize  lithographic 
patterning  at  the  nanoscale  to  transform  silicon  into  a 
material  with  an  effective  quantum  optical  nonlinearity  at 
an  engineerable  optical  wavelength.  The  modest  level  of 
squeezing  realized  in  this  work  is  predominantly  limited 
by  thermal  noise,  but  also  by  the  efficiency  with  which 
the  reflected  light  can  be  collected  by  external  optics. 
The  effects  of  thermal  noise  may  be  substantially  reduced 
by  working  with  materials  of  higher  intrinsic  Q- frequency 
product,  such  as  diamond  and  silicon  carbide  [29].  Given 
the  microchip  form  of  the  devices  studied  here,  and  the 
potential  for  device  integration,  it  is  interesting  to  con¬ 
sider  whether  new  squeezed  light  applications  might  arise. 
For  example,  squeezed  light  generated  by  one  device  could 
be  directly  sent  into  another  device  for  use  as  an  optical 
probe.  Such  an  on-chip  squeezer  and  detector  could  be 
used  as  a  quantum-enhanced  micro-mechanical  displace¬ 
ment  and  force  sensor  [30].  More  generally,  we  expect 
future  experiments  with  feedback  and  strong  measure¬ 
ment  of  the  dynamics  of  a  mechanical  system  to  be  within 
reach.  In  addition,  using  quantum  light  as  an  input  we 
expect  to  be  able  to  generate  entangled  states  of  mechan¬ 
ics  and  light  with  similar  devices. 
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I.  THEORY 

Optomechanical  systems  can  be  described  theoretically 
with  the  Hamiltonian  (see  main  text) 

H  =  fajj0a)a  +  hujmoWb  +  hgoa^aQ^  +  b),  (Al) 

where  a  and  b  are  the  annihilation  operators  for  pho¬ 
tons  and  phonons  in  the  system,  respectively.  Generally, 
the  system  is  driven  by  intense  laser  radiation  at  a  fre¬ 
quency  Co>l  5  making  it  convenient  to  work  in  an  interaction 
frame  where  uoQ  is  replaced  by  A  in  the  above  Hamilto¬ 
nian  with  A  =  uj0  —  co>l.  To  quantum  mechanically  de¬ 
scribe  the  dissipation  and  noise  from  the  environment,  we 


use  the  quantum-optical  Langevin  differential  equations 
(QLEs)  [26,  31,  32], 

a(t)  =  -  (iA  +  ^  a  -  igoafo  +  b) 

\/ ^edin  (t)  (^)  ? 

b(t)  =  -  (iwm0  +  b  -  ig0a)a  - 

which  account  for  coupling  to  the  bath  with  dissipation 
rates  /q,  Ke,  and  71  for  the  intrinsic  cavity  energy  decay 
rate,  optical  losses  to  the  waveguide  coupler,  and  total 
mechanical  losses,  respectively.  The  total  optical  losses 
are  k  =  Ke  +  /q.  These  loss  rates  are  necessarily  accom¬ 
panied  by  random  fluctuating  inputs  din(t),  «in,i  (t) ,  and 
bin(£),  for  optical  vacuum  noise  coming  from  the  coupler, 
optical  vacuum  noise  coming  from  other  optical  loss  chan¬ 
nels,  and  mechanical  noise  (including  thermal). 

The  study  of  squeezing  is  a  study  of  noise  propagation 
in  the  system  of  interest  and  as  such,  a  detailed  under¬ 
standing  of  the  noise  properties  is  required.  The  equa¬ 
tions  above  are  derived  by  making  certain  assumptions 
about  the  noise,  and  are  generally  true  for  the  case  of 
an  optical  cavity,  where  thermal  noise  is  not  present,  and 
where  we  are  interested  only  in  a  bandwidth  of  roughly 
108  smaller  than  the  optical  frequency  (0  -  40  MHz  band¬ 
width  of  a  200  THz  resonator).  For  the  mechanical  sys¬ 
tem,  where  we  operate  at  very  large  thermal  bath  oc¬ 
cupancies  (^>  103)  and  are  interested  in  the  broadband 
properties  of  noise  sources  (0  -  40  MHz  for  a  30  MHz 
resonator),  a  more  detailed  understanding  of  the  bath  is 
required,  and  will  be  presented  in  the  section  on  thermal 
noise. 

At  this  point,  we  linearize  the  equations  assuming  a 
strong  coherent  drive  field  cko,  and  displace  the  annihi¬ 
lation  operator  for  the  photons  by  making  the  transfor¬ 
mation  a  -A  <ao  +  a.  This  approximation,  which  neglects 
terms  of  order  a2  is  valid  for  systems  such  as  ours  where 
go  <C  /q  i.e.  the  vacuum  weak  coupling  regime.  We  are 
then  left  with  a  parametrically  enhanced  coupling  rate 
G  =  g0\a0\.  Using  the  relations  given  in  the  mathemati¬ 
cal  definitions  section  (VI)  of  this  document,  we  write  the 
solution  to  the  QLEs  in  the  Fourier  domain  as 

~  /  \  _  -q^inM  -  -  iG(b(u)  +  qq) 

“  *(A-q  +  «/2 

b(u)  =  _  iG(a(w)  +afM)  ^ 

i(wmo  -  q  +  7i/2  i(wm0  -  q  +  7i/2 ' 

Finally  we  note  that  by  manipulation  of  these  equa¬ 
tions,  the  mechanical  motion  can  be  expressed  as  a  (renor¬ 
malized)  response  to  the  environmental  noise  and  the  op¬ 
tical  vacuum  fluctuations  incident  on  the  optical  cavity 
through  the  optomechanical  coupling 
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b(u)  = 


i(wm  -  u)  +  7/2 

y7'<Ain  (^)  “1“  \/^'i^'in,i  (^) 

i(A  —  cj)  +  k/2  i(wm  —  uS)  +  7/2 

iG  V^«inH  +  V^Qjn.iH 

— i(A  +  c<;)  +  ft/2  i(c<;m  —  cj)  +  7/2 


(A3) 


The  renormalized  mechanical  frequency  and  loss  rate 
are  ujm  =  co’mo  #  &jm,  and  7  =  7i  +  70M,  respectively, 
with 


Suj  m  =  |G|2Im 
7om  =  2|G|2Re 


1 

^m)  T  ^/2 

1 

-  cjm)  T  ft/2 


1 

— i(A  +  cjm)  +  ft/2 
(A4) 

1 

— i(A  +  cjm)  +  ^/2 
(A5) 


It  is  convenient  to  define  here  what  we  mean  by  a 
quadrature,  as  it  is  the  observable  of  the  light  field  that 
our  measurement  device  (the  balanced  homodyne  detec¬ 
tor  (BHD)  setup)  is  sensitive  to: 

A'7  aje  i0,\  a)eie.  (A6) 

We  are  interested  in  the  properties  of  for  various 

quadrature  angles  #,  given  the  influence  of  the  mechanical 
system. 

The  measurement  of  the  field  provides  us  with  a  record 
f(t)  =  XgOUi\t )  for  a  certain  6.  We  use  a  spectrum  an¬ 
alyzer  to  perform  Fourier  analysis  on  this  signal  and  ob¬ 
tain  a  symmetrized  classical  power  spectral  density  (PSD) 
§n(uj ),  as  defined  in  the  mathematical  appendix  (section 
VI). 

For  a  vacuum  field  such  as  the  input  field,  the  measured 
quadrature  (t)  will  have  a  power  spectral  density 

7TM  =  1.  (A7) 


This  is  the  shot-noise  level  which  is  due  to  the  quan¬ 
tum  fluctuations  of  the  electromagnetic  field.  Mathemat¬ 
ically,  it  arises  from  the  correlator  (avac(co’)a^ac(co’/))  = 
S(uj  +  u/),  with  all  other  correlators  (d^ac(co’)dvac(co’/)), 
(4ac  7)^*7^  ’  (“vac  Mavac  (<*/)),  arising  in  the  expres¬ 
sion  (v{w)I(w'))  equal  to  zero. 


A.  Approximate  quasi-static  theory 

In  this  section  we  present  a  simplified  derivation  of  how 
squeezing  is  obtained  in  the  studied  optomechanical  sys¬ 
tem  to  elucidate  the  important  system  parameters  and 
their  role  in  squeezing.  We  make  a  few  approximations 
to  simplify  the  derivation: 


1.  A  =  0:  The  laser  is  tuned  exactly  to  the  optical 
cavity  frequency. 

2.  fte  =  ft:  Perfect  coupling. 

3.  ft  cjm:  Bad  cavity  limit. 

4.  uj  <C  cjm:  We  are  only  interested  in  the  quasi- static 
response,  so  the  resonant  response  of  the  mechanical 
resonator  does  not  play  a  role. 

Under  these  assumptions,  equations  (2)  and  (3)  can  be 
written  as  (using  the  relation  for  the  optical  output  field 

flout  M-  =  ain(uj)  +  v^aM): 

zwm6(w)  =  -VTi^in  (w)  +  ^k(ain(w) +a|„(w)), 

V  K 

flout  M'  =  -flm(w)  -  +  6f(w)).  (A8) 

v  K* 

The  first  equation  shows  the  response  of  the  mechani¬ 
cal  resonator  subsystem  to  the  thermal  bath  fluctuations 
(b[n(uj))  and  the  optical  vacuum  noise  from  the  measure¬ 
ment  back-action.  We  define  rmeas  =  4|G|2/ft,  and  inter¬ 
pret  it  as  the  measurement  rate  [2],  such  that  the  factor 
appearing  in  front  of  the  optical  vacuum  noise  operators 
is  V/Tmeas-  This  rate  also  appears  in  the  second  equation 
for  the  output  field,  in  front  of  the  normalized  position 
operator  x/xzp{  =  b(oj)  +  which  is  the  observable 

that  is  being  measured. 

Note,  from  the  expression  for  aou t(cj)  it  follows,  that 
since  the  position  is  a  real  observable  with  an  imaginary 
prefactor,  the  effects  we  consider  depend  strongly  on  the 
quadrature  being  probed,  i.e.  the  real  part  of  the  expres¬ 
sion,  Xq°^ ,  will  not  be  affected  by  the  optomechanical 
coupling. 

At  this  point  we  can  easily  calculate  the  properties  of 
the  detected  spectrum  Sjft( u),  by  writing  dout  in  terms 
of  din  and  b[n  for  which  the  correlators  are  known: 

flout (w)  =  -OinM  ~  2i'^meas  (ain(a?)  +  ajnM) 

+  -»!>)).  (AS!) 

idm 

Ignoring  thermal  noise  for  the  moment  (71  =  0),  and  drop- 
ping  terms  of  order  (rmeas/wm)2  (assuming  rmeas  <  wm) 
we  arrive  at: 


=  1 +  4(rmeas/wm)sin(26>).  (A10) 

Note  that  for  certain  values  of  #,  the  detected  spectral 
density  can  be  smaller  than  what  one  would  expect  for  a 
vacuum  field.  For  6  =  — 7r/ 4,  we  achieve  the  maximum 
squeezing  with  a  noise  floor  of  1  —  4(rmeas/cjm)  which 
strongly  dependends  on  the  ratio  Tmeas/um. 

To  understand  the  effect  of  thermal  noise,  we  as¬ 
sume  the  form  of  the  correlator  to  be  (6in  M^„<y)>  = 
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Figure  Al:  Squeezing  theory,  a,  Density  plot  of  the  predicted  squeezing  *§//t(cu)  vs.  phase  angle  and  frequency,  normalized 
to  the  shot-noise.  The  mechanical  mode  can  clearly  be  seen  at  cum/27r  =  30  MHz.  The  solid  white  lines  outline  the  region 
where  the  power  spectral  density  falls  below  1  (the  shot-noise  level)  indicating  the  presence  of  squeezing  for  that  phase  and 
frequency.  The  dashed  white  lines  at  0  =  — 7t/4  and  6  —  +7t/4  correspond  to  regions  where  squeezing  can  be  obtained  below 
and  above  the  mechanical  frequency,  respectively,  and  the  components  of  the  noise  model  for  these  phases  is  shown  in  detail 
in  figures  b  and  c.  In  these  figures  the  spectra  are  again  normalized  to  the  shot-noise  level  plotted  as  a  grey  line.  The  simple 
squeezing  model  without  thermal  noise  (Eq.  (A10))  is  represented  by  the  dashed  green  line  and  the  simple  model  with  thermal 
noise  (Eq.  (All))  is  the  solid  green  line.  The  solid  black  line  is  the  full  squeezing  model  5//t(cu)  corresponding  to  a  with 
the  constituent  components:  the  contribution  from  the  optical  vacuum  fluctuations  (Sj^a(oj);  Eq.  (A24))  represented  by  the 
dashed  black  line  and  the  thermal  noise  Eq.  (A25))  represented  by  the  dashed  red  line. 


(jn{uS)  +  1)5(cj  +  cu'),  (bln(id)b[n(uj'))  —  n(uj)S(u:  +  (*/), 
(bln(u:)bln(uj'))  =  0,  and  (bin(w)bin(w'))  =  0  (these  expres¬ 
sions  are  discussed  in  section  IV  G).  Then  a  calculation 
similar  to  the  one  leading  to  equation  (10)  gives 


main  text,  elucidates  the  role  of  correlations  between  the 
mechanical  system’s  position  and  the  back-action  force. 

The  response  of  the  mechanical  system  to  a  force  is 
captured  by  its  susceptibility: 


=  1  +  4(rmeasMn)sin(20) 

+4rmea8y  (1  +  cos(2g))>  (A11) 

where  we  have  assumed  n(cj),  the  bath  occupation  at  fre¬ 
quency  cj,  to  be  much  larger  than  unity.  At  0  =  — 7r / 4, 
we  have 

=  1  -  4(rmeasMn)(l  -  n(u)/Qm).  (A12) 

In  this  model,  there  is  no  squeezing  at  c o  if  n(oo)  >  Qm. 
This  means  that  for  squeezing  to  be  present,  coherent 
evolution  of  the  mechanical  system  must  be  possible,  i.e. 
the  rate  at  which  phonons  enter  the  mechanical  system 
from  the  bath  (71  n)  must  be  smaller  than  the  mechanical 
frequency  wm.  In  conclusion,  the  important  requirements 
to  achieve  squeezing  are  to  make  rmeas  comparable  to 
cjm,  and  to  reduce  the  thermal  occupancy  or  increase  the 
mechanical  quality  factor  to  achieve  n{uS)  <  Qm. 


B.  The  effect  of  dynamics  and  correlation  between 
RPSN  and  position 


XmM 


1 

m(u4  -  u!2  -  i7iWm)  ’ 


(A13) 


The  form  of  the  damping  considered  here  is  the  strongly 
sub-ohmic  structural  damping  which  is  observed  in  our 
measurements  [33,  34]  (cf.  Section  IV  G).  The  mechanical 
system  responds  to  random  noise  forces  FT(t)  from  the 
thermal  bath  (which  we  treated  in  the  last  section  and 
neglect  here),  and  to  the  quantum  back-action  from  the 
cavity  FBa(*)- 

The  back-action  force  for  the  resonant  case  can  simply 
be  found  by  linearizing  the  expression  for  the  radiation 
pressure  force  Frp(£)  =  —hga^a/xz pf.  We  find  the  force 
imparted  on  the  mechanics  due  to  the  shot-noise  of  the 
cavity  field  to  be 


PBA(t)  =  (A14) 

^zpf 

for  the  case  of  resonant  driving.  The  fluctuations  im¬ 
parted  on  the  mechanics  are  from  the  intensity  quadra¬ 
ture  of  the  light  (6  =  0).  Using  equation  (8),  we  can  write 
the  output  field  quadrature  as: 


As  a  next  step,  we  take  into  account  the  dynamics  of  the 
mechanical  resonator  while  keeping  the  approximations  of 
the  bad-cavity  limit  (k  oom)  and  on-resonant  probing 
(A  =  0).  In  addition  to  further  clarifying  some  of  the 
observed  features,  this  treatement,  as  presented  in  the 


xfut\t)  =  -x£n\t)  -  2 -  sin(0).  (A15) 

^zpf 

We  note  here  that  the  mechanical  position  fluctuations 
are  primarily  imprinted  on  the  phase  quadrature  of  the 
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output  light,  with  0  =  ±tt/2.  The  intensity  quadrature 
is  unmodified  =  —  X^0{t))  since  changes  in 

the  cavity  frequency  are  not  transduced  as  changes  in 
intensity  when  the  laser  is  resonant  with  the  cavity. 

The  output  of  the  homodyne  detector  normalized  to  the 


shot-noise  level  is  found  by  taking  the  auto-correlation  of 
eqn.  (A15).  The  correlations  between  radiation  pressure 
shot-noise  and  the  mechanical  motion  are  important  in 
this  calculation  [35-41]  and  must  be  taken  into  account. 
In  the  time-domain  we  find  the  auto-correlation  to  be: 


<x(out)(^OUtV)>  =  ^-0+4rmeas  sin2(0)^kM 

Xzpf 

+2 hr1  sin(6^)  cos(6)(FBA(t)x(t')  +  (A16) 


The  cos($)  in  the  last  term  comes  from  the  general  ex¬ 
pression  for  a  quadrature  X^n\t)  =  X^?q(£)cos($)  + 
X^\2(t)  sin($),  and  equation  (A14).  The  key  compo¬ 
nents  of  equation  (A16)  are  the  shot-noise  level,  the  ther¬ 
mal  noise,  and  the  cross-correlation  between  the  back- 
action  noise  force  and  mechanical  position  fluctuations. 
It  is  only  the  latter  which  can  give  rise  to  squeezing, 
by  reducing  the  fluctuation  level  below  shot-noise.  This 
squeezing  can  be  calculated  spectrally: 

/oo 

dr[(FBA(i)^(i  -  t)) 

-OO 

+(x(t)FBA(t  -  T))]elbJT 

=  2^sin(26»)l^^xm(w) 

^zpf 

=  4mwmsin(20)rmeasXm(w).  (A17) 

At  the  DC  or  quasi-static  limit  (cu  0)  the  susceptibility 
Xm  — >*  Vmc°m  can  be  used  and  we  reobtain  the  results 
from  section  I A  (cf.  equation  (A10)).  We  see  that  for 
0  <  0,  squeezing  is  obtained  in  this  limit.  At  frequencies 
larger  than  u;m,  Xm(^)  changes  sign,  and  we  expect  to  see 
squeezing  at  quadrature  angles  0  >  0.  Additionally,  since 
Xm(w)  becomes  larger  around  the  mechanical  frequency, 
we  expect  the  maximum  squeezing  to  be  enhanced.  More 
specifically,  at  a  detuning  S  =  oo  from  the  mechanical 
resonance,  we  expect  the  parameter  characterizing  the 
squeezing  to  be  proportional  to  Tmeas/(k  Ah  of  these 
features  are  seen  in  Fig.  1. 

It  is  important  to  note  here  that  in  the  absence  of  other 
nonlinearities  in  the  system,  any  reduction  of  the  noise 
below  the  vacuum  fluctuations  can  only  be  caused  by  the 
correlations  between  the  RPSN  and  the  position  fluctua¬ 
tions  of  the  system.  This  makes  the  problem  of  proving 
the  correlations  between  RPSN  and  mechanical  motion 
equivalent  to  the  problem  of  proving  that  the  reflected 


light  from  the  optomechanical  cavity  has  been  squeezed. 

Conceptually  this  form  of  probing  the  RPSN  is  similar 
to  that  carried  out  by  Safavi-Naeini  et  al.  [38,  41]  and  ana¬ 
lyzed  by  Khalili  et  al.  [39].  It  also  shares  features  with  the 
cross-correlation  measurements  proposed  by  Heidmann  et 
al.  [35],  and  Bprkje  et  al.  [37],  and  recent  experiments  by 
Purdy  et  al.  [6].  The  distinguishing  feature  of  this  type 
of  measurement  is  that  the  quantum  correlations  between 
the  fluctuations  of  the  position  and  the  electromagnetic 
vacuum  manifest  themselves  as  squeezed  light. 


C.  General  derivation  of  squeezing 

Among  the  approximations  made  in  section  I  A,  the 
quasi- static  approximation  is  the  least  correct.  In  fact, 
in  our  experiments,  the  most  observable  squeezing  occurs 
with  l jj  close  to  ojm  and  even  slightly  larger  than  co>m,  so 
oo  <C  oo m  is  not  valid.  Near  the  mechanical  frequency, 
resonant  enhancement  of  the  optical  vacuum  fluctuations 
by  the  mechanical  resonator  causes  squeezing  greater  than 
that  predicted  in  the  quasi-static  regime  to  be  possible. 

Here  we  show  the  results  of  a  derivation  that  does  not 
rely  on  most  of  the  assumptions  used  in  the  approximate 
model.  Of  the  assumptions  in  the  previous  section,  the 
only  simplification  we  keep  here  is  to  assume  perfect  cou¬ 
pling  Ke  =  k.  The  effect  of  imperfect  coupling  can  be 
taken  into  account  trivially  and  is  explained  after  this 
section  (see  IC  1). 

By  substitution  of  equation  (3)  into  the  equation  for 
a(cu)  (2),  we  arrive  at: 

\Zk&(uj )  =  Ai(w)ain(w)  +  A2(u)a\n(w) 

+.Bi(w)Sjn(w)  -|-  i?2(w)S|n(a;),  (A18) 

with 
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Ai(u) 


A2(uj) 


B^lj) 

B2{lS) 


K 

i(  A  —  uS)  +  k/2 


|G|2 

1 

i( A  —  uj)  +  k/2  i(u) m  —  oo)  +  7/2 

K 

i(  A  —  uS)  +  k/2 

z\ 

to 

1 

[— i(A  +  uf)  +  k/2  i(cjm  —  uj)  +  7/2 

iG 

i{  A  —  w)  +  k/2 

_i(wm  -  cj)  +7/2_ 

iG 

i(  A  —  uS)  +  k/2 

_—i(cdm  +  cu)  +  7/2 

\G\2 _ 1  ] 

i(  A  —  cj)  +  k/2  —i(uj  m  +  cj)  +  7/2 

|G|2 _ 1 

— i(A  +  c j)  +  k/2  —  i(cjm  +  cj)  +  7/2 


(A19) 


(A20) 

(A21) 

(A22) 


The  expressions  give  us  the  output  field  in  terms  of  the 
input  fields,  since 


^out(^)  —  din(cj)  T  v^M) 

=  (1  +  Ai(u))ain(w)  +  A2(cj)aJn(cc;) 

+^l(+)frin(+)  +  -^2  (^)  &in  (^)  •  (A23) 


We  can  calculate  from  this  expression,  which  we 

split  into  two  parts,  one  only  due  to  the  optical  vacuum 
fluctuations,  and  the  other  containing  the  contribution 
from  thermal  noise:  =  Sjf*a(uj)  +  S/^b(u). 


^/qat+O  ~  |A2(— cj)|2  +  |1  +  Ai(u;)|2  +  2Re{e  2ie[\  +  Ai(lS))A2(— cj)}  (A24) 

SZtM  =  |Si(^)|2(fi(a7)  +  1)  + 

+  \B2(-u)\2(n(u)  +  1)4  \B2(u)\2n(u) 

+2Re{e-2i^i(cj)52(-o;)}(n(cj)  +  1)  +  2Re{e~2ie  B^-^B^jn^)  (A25) 


1.  The  effect  of  imperfect  optical  coupling  and  inefficient 
detection 

At  every  juncture  in  an  experiment  where  the  opti¬ 
cal  transmission  efficiency  is  less  than  unity  (77  <  1), 
an  equivalent  optical  circuit  can  be  defined  involving  an 
r]  :  (1  —  7)  beam  splitter  with  the  output  being  7  times  the 
input  and  (1  —  77)  times  the  vacuum.  Therefore  the  effect 
of  optical  losses  and  coupling  inefficiencies  on  the  detected 
spectra  can  be  calculated  by  replacing  the  measured  field 
quadrature  with: 

+det)  =  y^+out)  +  y/l  -  ??+vac)  (A26) 

This  source  of  vacuum  noise  is  completely  unrelated  to  the 
cavity  output,  and  there  are  no  cross-correlation  terms, 
so  the  detected  current  spectral  density  will  be  given  by 

7f  M  =  VSP\u)  +  (1  -  V)Sff(w),  (A27) 

where  SJjc(uj)  =  1  is  the  shot-noise. 

Measurement  inefficiencies  take  two  forms,  one  is  due  to 
ineffeciencies  in  the  detection,  while  the  second  is  because 


of  excess  electronic  noise  or  “dark  current”  present  due 
to  the  circuitry  of  the  detector  and  amplifier.  This  excess 
noise  can  also  be  thought  of  as  a  detection  inefficiency 
by  considering  the  amount  of  optical  shot-noise  inserted 
into  the  signal  which  would  produce  it.  Since  the  dark- 
current  is  measured  with  no  optical  input,  and  the  real 
shot-noise  level  increases  linearly  with  the  local  oscillator 
(LO)  power,  this  inefficiency  is  power  dependent  and  can 
be  minimized  for  large  LO  powers.  In  our  case,  the  dark 
current  was  found  to  be  10.4  dB  below  the  detected  shot- 
noise.  The  total  detector  efficiency  was  measured  to  be 
7hd  —  66%. 


II.  EXPERIMENT 

A.  Measurement  of  losses 

In  order  to  estimate  the  total  squeezing  expected  in 
our  setup  we  carefully  characterize  all  losses  in  our  sys¬ 
tem.  Some  of  these  losses  are  static  (e.g.  circulator  losses) 
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Figure  A2:  Experimental  setup.  A  detailed  description  of  the  experimental  setup  can  be  found  in  the  main  text. 


while  others  can  vary  from  experiment  to  experiment  (e.g. 
coupling  efficiency  of  the  fiber  taper  to  the  waveguide). 
In  figure  2  typical  losses  are  shown  as  efficiencies  (77)  for 
various  parts  of  the  experiment.  The  efficiency  of  send¬ 
ing  light  from  port  1  to  2  of  our  optical  circulator  is 
7712  =  85%,  and  7723  =  88%  for  port  2  to  3.  In  addition, 
the  efficiency  from  port  3  of  the  circulator  to  the  homo¬ 
dyne  detector  is  773H  =  92%.  All  these  losses  are  fixed  and 
do  not  change  over  time  as  the  components  are  optically 
spliced  together.  Measuring  the  coupling  efficiency  of  the 
fiber  taper  to  the  waveguide  is  done  every  time  a  new  data 
set  is  taken.  This  is  accomplished  by  switching  the  light 
that  is  reflected  from  the  waveguide  to  a  power  meter  and 
comparing  the  reflected  power  to  a  known  input  power 
with  the  laser  tuned  off-resonance  from  the  optical  mode 
(off-resonance  the  device  acts  as  a  near-perfect  mirror). 
Typical  achieved  efficiencies  are  around  77CP  =  90%.  The 
efficiency  of  the  homodyne  detection  strongly  depends  on 
the  alignment  of  the  polarization  between  the  local  oscil¬ 
lator  and  the  signal,  as  well  as  by  how  much  the  power  in 
the  LO  overcomes  the  electronic  noise  floor  of  the  detec¬ 
tor.  To  determine  this  efficiency  we  use  an  acousto-optic 
modulator  (AOM)  inserted  in  our  setup  before  the  circu¬ 
lator  in  the  signal  path.  The  AOM  shifts  the  frequency 
of  the  light  creating  a  tone  88  MHz  away  from  the  signal 
with  a  fixed,  known  amplitude,  and  identical  polariza¬ 
tion  to  the  signal  (we  directly  measure  the  power  of  this 
tone  with  a  power  meter).  This  tone  can  now  be  used 
to  determine  the  total  homodyne  efficiency  by  measur¬ 
ing  its  power  on  the  spectrum  analyzer,  taking  the  other 
losses  into  account.  Our  typical  homodyne  efficiency  is 
77hd  =  66%  resulting  in  a  total  setup  efficiency  (detection 
efficiency  of  optical  signal  photons  in  the  on-chip  waveg- 
uide)  of  roughly  rjSe tup  =  Vgp  •  V23  ■  %h  •  hHD  ~  48%. 


B.  Data  collection  procedure 

Careful  calibration  of  our  data  is  crucial  in  understand¬ 
ing  all  noise  sources  and  potential  drifts  over  time  in  our 
setup.  The  losses  in  our  setup  are  determined  before  we 
make  a  new  data  run  as  described  in  the  previous  section. 


0  0.01  0.02  0.03  0.04  0.05 

Time  (s) 

Figure  A3:  Phase  information.  The  blue  trace  shows  the 
interference  signal  of  the  local  oscillator  and  the  signal  on  the 
homodyne  detector  when  their  relative  phase  is  scanned  using 
a  fiber  stretcher  in  the  local  oscillator  (LO)  arm.  The  voltage 
reading  here  is  proportional  to  cos(0  —  (f>)  where  0  is  the  phase 
difference  between  the  LO  and  input  to  the  cavity,  and  0  is 
the  phase  imparted  by  reflection  off  the  cavity.  This  inter¬ 
ference  signal  is  used  to  actively  stabilize  the  relative  phase 
to  different  set  points  (green  traces).  Occasionally  the  lock 
fails,  as  shown  by  the  red  traces,  and  any  associated  data  is 
discarded.  The  range  in  which  the  phase  can  be  stably  locked 
is  slightly  smaller  than  —  n/2  to  n/2  due  to  the  turning  points 
in  the  sinusoidal  interference  curve. 


We  then  proceed  to  record  an  optical  trace  of  the  cavity 
resonance  by  switching  the  light  to  a  photodetector  (PD1 
in  figure  2)  and  scanning  the  laser  wavelength.  This  trace 
provides  the  information  to  lock  the  laser  to  a  fixed  detun¬ 
ing  (typically  0.04  •  k  red  of  the  cavity  resonance),  which 
is  accomplished  using  a  simple  software  lock  and  feed¬ 
back  from  the  wavemeter  (with  a  resolution  of  roughly 
0.003  •  k)  and  is  described  in  more  detail  in  the  subsec¬ 
tion  below.  As  a  next  step  the  optical  signal  is  switched 
to  the  homodyne  detector  and  the  relative  phase  between 
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the  signal  and  the  local  oscillator  is  scanned  using  the 
fiber  stretcher  in  the  LO  arm.  The  resulting  interference 
is  shown  in  figure  3  as  the  blue  trace.  The  interference 
signal  is  used  to  lock  the  relative  phase  between  the  signal 
and  LO  using  a  Toptica  DigiLock  110.  The  green  traces 
show  the  properly  locked  signal,  while  the  red  traces  are 
phase  set  points  where  the  lock  failed  requiring  the  asso¬ 
ciated  data  to  be  discarded.  We  then  record  the  spectra 
of  the  homodyne  signal  and  for  every  trace  taken  we  also 
save  a  spectrum  of  the  shot-noise  by  switching  the  signal 
arm  away  from  the  homodyne  detector  and  only  measure 
vacuum  input  to  the  signal  arm  of  our  detector.  We  re¬ 
lock  the  laser  with  respect  to  the  cavity  every  other  data 
point  to  counteract  drift.  This  procedure  is  repeated  for 
several  different  phases  and  different  input  powers.  We 
typically  took  data  for  60  different  phases  for  every  input 
power  within  a  range  of  a  little  less  than  — 7r/2  to  7r/2. 


C.  Relation  between  detuning  and  quadrature 


The  laser  frequency  is  positioned  at  a  detuning  of 
roughly  0.04-  k  by  starting  at  a  larger  detuning  on  the  red 
side  of  the  cavity,  and  stepping  the  laser  blue  in  0.1  pm 
steps  (12  MHz)  towards  the  cavity  while  monitoring  the 
average  intensity  of  the  reflected  light  on  PD1.  Once 
the  target  intensity  is  reached,  the  laser  is  kept  at  this 
wavelength  during  the  course  of  the  measurement  by  the 
wavemeter  lock  without  further  feedback  from  PD1.  The 
intensity  reading  gives  us  an  idea  of  the  value  of  the  de¬ 
tuning  which  is  determined  more  accurately  by  analysis 
of  our  homodyne  spectra. 

The  homodyne  spectra  are  taken  at  different  phase  lock 
points  (see  Fig.  3)  corresponding  to  quadrature  angles 
0\ock  between  the  reflected  signal  and  the  LO.  These  an¬ 
gles  differ  from  our  convention  in  Section  I  where  the 
phase  0  between  the  input  light  into  the  cavity  and  the 
local  oscillator  is  considered.  They  are  related  to  onean- 
other  by  the  phase  imparted  on  the  input  light  upon  re¬ 
flection  from  the  cavity, 


4>(A)  =  Arg 


HZe 

zA  T  k/2 


and  the  relation 


(A28) 


tfiock  =  0-0.  (A29) 

For  a  given  laser-cavity  detuning  A,  we  sweep  through 
the  different  phase  lock  points  (see  Fig.  3)  #iock,  and 
take  mechanical  spectra  for  each  phase.  The  phase  that 
minimizes  the  mechanical  signal  9*ock  is  determined  from 
the  recorded  spectra.  This  allows  us  to  solve  for  A  us¬ 
ing  the  expression  0*ock  =  0*( A)  —  0(A)  where  0*( A) 
is  the  phase  minimizing  the  mechanical  transduction  ac¬ 
cording  to  the  model  in  the  previous  section.  To  first 
order  (for  A  <C  k)  0*  is  0  since  no  mechanical  signal 
is  observed  in  the  intensity  quadrature  of  the  reflected 
light.  This  post-processing  of  the  data  allows  us  to  deter¬ 
mine  that  across  the  measured  powers  the  detuning  was 


Normalized  detuning,  A/k 


Figure  A4:  Detuning  and  phase  lock  points,  a,  An  optical 
scan  taken  before  the  data  run  starts  is  shown.  The  blue  ver¬ 
tical  line  denotes  the  target  detuning  the  software  lock  moves 
the  laser  to,  determined  from  the  measured  reflection  inten¬ 
sity.  The  laser  is  kept  at  that  detuning  via  a  wavemeter  lock, 
as  the  light  is  switched  away  from  PD1  and  to  the  homodyne 
detector.  The  measured  area  under  the  mechanical  mode  is 
plotted  in  b  (blue  circles)  at  this  detuning.  A  minimum  value 
is  reached  for  a  local-oscillator  to  reflected  signal  phase  of  0*ock. 
Depending  on  the  detuning,  different  mechanical  mode  ampli¬ 
tudes  can  be  measured  at  this  phase  angle  #iock,  according  to 
the  model.  We  obtain  an  accurate  estimate  of  the  detuning  by 
calculating  the  detuning  at  which  the  mechanical  mode  am¬ 
plitude  is  minimized  at  the  measured  0\Ock  as  shown  in  c.  The 
expected  mode  amplitudes  for  the  detunings  represented  by 
the  red  and  green  lines  in  c  are  shown  by  the  similarly  colored 
curves  in  b. 


A  =  0.044  •  k  d=  0.006  •  k.  For  a  single  measured  power, 
we  expect  a  more  accurate  determination,  with  an  uncer¬ 
tainty  on  the  order  of  0.003  •  k.  This  level  of  accuracy  in 
the  detuning  also  determines  the  uncertainty  in  quadra¬ 
ture  angle  of  0.04  rad. 


III.  SAMPLE  FABRICATION  AND 
CHARACTERIZATION 


A.  Fabrication 

The  devices  are  fabricated  from  a  silicon  on  insulator 
(SOI)  wafer  (SOITEC,  220  nm  device  layer,  3  fim  buried 
oxide,  device  layer  resistivity  4  —  20  £4  •  cm)  using  elec¬ 
tron  beam  lithography  and  subsequent  reactive  ion  etch¬ 
ing  (RIE/ICP)  to  form  the  structures.  The  buried  oxide 
is  then  removed  in  hydrofluoric  acid  (49%  aqueous  HF  so¬ 
lution)  and  the  devices  are  cleaned  in  a  piranha  solution 
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(3:1  H2SO4  and  H2O2)  and  finally  hydrogen  terminated 
in  diluted  HF.  For  a  more  detailed  description  see  [42]. 

B.  Optical  Characterization 

The  optical  characterization  of  our  devices  is  done 
by  sweeping  the  laser  frequency  across  the  optical  res¬ 
onance  while  detecting  the  reflected  light  in  a  photode¬ 
tector  (PD1  in  figure  2).  This  light  is  simultaneously  sent 
to  a  wavemeter  to  record  the  absolute  wavelength  and  ac¬ 
curately  determine  the  linewidth  and  center  frequency  of 
the  resonance.  Each  chip  contained  several  designs  where 
the  waveguide  loading  (coupling)  of  the  optical  cavity  was 
varied  by  changing  the  gap  size  between  the  waveguide 
and  nanobeam.  For  our  measurements  we  chose  a  slightly 
overcoupled  (fte//q  ~  1.22  >  1)  device  with  good  optical 
quality  (57,000  loaded  Q)  [43]. 

C.  Mechanical  Characterization 


The  intrinsic  mechanical  damping  rate  71  and  the  op¬ 
tomechanical  coupling  rate  go  are  measured  by  detecting 
the  mechanical  response  to  the  signal  laser,  through  the 
reflected  signal  field,  on  the  spectrum  analyzer.  We  keep 
the  optical  power  constant,  while  we  take  measurements 
for  several  different  detunings  A.  The  radiation  pres¬ 
sure  force  causes  both  an  optical  spring  effect  resulting 
in  a  frequency  shift  of  the  mechanical  resonance,  as  well 
as  damping  of  the  mechanical  motion,  associated  with  a 
broadening  of  its  linewidth  (see  equations  (A4)  and  (A5)). 
By  fitting  the  data  shown  in  figures  5a  and  5b,  we  can 
extract  71  =  172  Hz  and  go  =  27r  x  750  kHz.  Knowing 
the  mechanical  properties  of  our  resonator  and  the  pre¬ 
cise  intracavity  photon  number,  we  can  now  also  extract 
the  mechanical  bath  occupancy  no  as  a  function  of  de¬ 
tuning  from  the  measured  PSDs  of  the  mechanical  mode 
(figure  5c;  see  also  [24]).  This  shows  us  that  our  mechan¬ 
ical  mode  thermalizes  to  about  16  K  for  low  optical  input 
powers,  which  is  close  to  the  cold  finger  temperature  of 
our  cryostat  of  10  K. 

IV.  NOISE  SPECTROSCOPY  DETAILS 
A.  Homodyne  measurement  with  laser  noise 

Our  experiment  is  designed  to  measure  the  spectral 
density  of  the  fluctuations  of  the  optical  field  exiting  the 
cavity.  However  any  real  laser  system  will  have  technical 
noise,  in  addition  to  the  quantum  noise  associated  with 
an  ideal  coherent  source,  which  adds  to  the  detected  noise 
level.  Both  the  signal  and  local  oscillator  arm  of  our  setup 
contain  this  noise  which  must  be  taken  into  account.  The 
noise  on  the  signal  arm  can  also  be  modified  non-trivially 
by  propagation  through  the  optomechanical  system.  We 


start  by  reproducing  known  results  on  the  operation  of  an 
ideal,  balanced  homodyne  detection  system  with  signal 
and  local  oscillator  input  fields  as  and  a^o  respectively, 
under  the  influence  of  noise  [28,  44,  45].  Most  generally, 
these  fields  consist  of  coherent  tones  as  and  070  >  techni¬ 
cal  (or  classical)  noise  components  as^(t)  and  aoo^(t)} 
and  quantum  fluctuations  as?vac(t)  and  nLO,vac(^): 

CLS  =  Cys  T  Us,n(^)  "b  ^s,vac(^)}  (A30) 

«LO  =  OLO  +  «lO,nW  +  dLO,vac(£)-  (A31) 

Since  both  the  local  oscillator  field  and  the  signal  field 
are  generated  by  the  same  laser,  the  technical  noise  on 
the  signal  and  local  oscillator  will  be  correlated,  and  these 
correlations  must  be  accounted  for  in  the  analysis.  For  the 
simplest  case,  where  the  signal  arm  does  not  experience 
the  complex  dispersion  from  interaction  with  an  optome¬ 
chanical  system  (e.g.  being  reflected  off  the  end-mirror  far 
detuned  from  the  optical  resonator),  we  expect 

as,N(i)  =  aB£(t)  and  aLO,N  (t)  ~  aLO£(t).  (A32) 

The  function  £(t)  is  related  to  the  intensity  and  phase 
fluctuations  of  the  laser  light  (n(t)  and  </>(£)  respectively): 

a(t)  =  a0(l  +  n(t))e1^  &  a0(l  +  n(t)  +  i(j>(t)) 

£(t)  =  n(t)  +  i(j)(t)  (A33) 

The  difference  of  the  photocurrent  in  the  homodyne 
detector  is  given  by 

I  (/)  =  as®LO  +  «LO ,  ( A34) 

which,  considering  only  the  technical  noise,  reduces  to 

I{t)  =  |aLO|^s’vac)  +  /dc(1  +  2Re{£(*)}),  (A35) 

under  the  assumption  that  a^o  using  the  def¬ 

initions  in  equation  (32),  and  taking  the  DC  current 
Jdc  =  2Re{<a*aLo}  =  2|<asa'Lo|  cos(0),  where  6  is  the  rel¬ 
ative  phase  between  the  signal  and  local  oscillator.  From 
this  equation  we  see  that  the  phase  noise  <j>(t)  cannot  be 
detected  on  a  balanced  homodyne  setup.  This  can  be 
understood  as  being  from  the  detectors  fundamental  in¬ 
sensitivity  to  phase  noise  on  the  laser,  as  the  only  phase 
reference  in  the  system  is  the  local  oscillator,  which  con¬ 
tains  the  same  phase  fluctuations  as  the  signal.  Secondly, 
for  the  local  oscillator  phase  which  makes  Jdc  =  0,  inten¬ 
sity  noise  is  not  detected.  In  a  real  homodyne  detector 
this  is  only  true  for  a  perfect  common  mode  rejection 
ratio  (CMRR),  which  is  the  case  in  our  setup  as  the  in¬ 
tensity  noise  is  negligible  and  the  CMRR  is  >  25  dB. 
For  these  reasons  we  use  a  different  setup  for  character¬ 
izing  the  laser  phase  and  intensity  noise  as  described  in 
Section  IV  B. 

B.  Measurement  and  characterization  of  laser  noise 

In  this  section  we  discuss  the  procedure  used  for  charac¬ 
terization  of  our  laser  (New  Focus  TLB-6728-P-D).  This 
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Figure  A5:  Optomechanical  characterization.  We  characterize  the  behavior  of  the  optomechanical  system  in  order  to 
extract  several  parameters  such  as  the  intrinsic  mechanical  linewidth  71,  the  optomechanical  coupling  rate  go ,  and  the  bath 
temperature  To  (n b).  a,  The  effective  mechanical  frequency  =  a;mo  +  8w  described  in  equation  (A4)  is  plotted  as  a  function 
of  the  laser  detuning  A  =  uj0  —  ujo  (shown  here  in  units  of  wavelength) .  The  frequency  shift  is  due  to  the  optical  spring  effect 
caused  by  radiation  pressure,  b,  The  optomechanical  interaction  also  causes  the  intrinsic  linewidth  71  of  the  mechanical  mode 
to  be  broadened  as  the  detuning  is  changed  (cf.  equation  (A5)).  c,  The  area  under  the  mechanical  Lorentzian  is  also  modified 
depending  on  A,  and  is  shown  here,  normalized  to  shot-noise.  The  fits  (green  lines)  in  a-c  are  now  used  to  obtain  71,  go  and 
rib  (see  text  for  details).  The  plot  in  d  shows  a  normalized  cavity  scan,  which  is  used  to  determine  the  exact  detunings  in  a-c, 
with  every  red  data  point  corresponding  to  a  data  point  in  a,  b  and  c. 


For  the  phase  noise 
measurement 


Figure  A6:  Experimental  setup  for  characterization  of 
intensity  and  phase  noise.  The  laser  is  amplitude  stabilized 
and  an  attenuator  is  used  to  select  the  desired  optical  power. 
For  the  phase  noise  measurement  the  light  is  sent  through  a 
Mach-Zehnder  interferometer  (MZI)  with  a  free  spectral  range 
of  115  MHz.  The  laser  is  locked  to  the  center  of  the  interfer¬ 
ence  fringe  allowing  frequency  noise  to  be  converted  to  inten¬ 
sity  noise.  The  light  is  then  detected  on  a  New  Focus  Model 
1811  photodetector  and  the  photocurrent  detected  on  a  spec¬ 
trum  analyzer.  The  same  setup  is  used  to  detect  intensity 
noise  without  the  MZI. 


characterization  was  done  using  an  independent  setup, 
shown  in  figure  6,  and  involved  two  measurements  di¬ 
rectly  detecting  the  light. 

The  first  measurement  is  to  characterize  the  intensity 
noise  where  the  laser  light  is  sent  directly  onto  a  photode¬ 
tector  with  the  incident  power  varied.  From  the  theory 
we  expect  for  the  detector  photocurrent 

I(t)  =  (<aLo  +  Rlo,n(£)  +  ^lo  ,vac  W)t(  h.c.) 

=  |aLo|y£°’vac)  +  W1  +  2R e{£(i)}), 
with  JDC  =  |aLo|2- 

The  spectral  density  of  the  current  is  then  given  by 

Sll(id)  =  |c^LO  |2  (l  +  |«Lo|  (gj))  ,  (A36) 


where  Snn(uj )  is  the  PSD  of  the  intensity  noise  fluctua¬ 
tions.  For  a  real  detector,  this  equation  is  modified  by 
the  presence  of  a  dark  current  Sdark(^)  and  non- unity 
efficiency  (igdet  <  1): 

Sn(u)  =  S'dark(^) 

+  |«Lo|2  (l  +  ^det|0T0|2  (gj))  .  (A37) 

We  subtract  the  dark  current  (measured  with  the  laser 
turned  off)  from  the  data,  and  set  bounds  on  the  mag¬ 
nitude  of  the  intensity  noise  present  in  the  laser  by  ex¬ 
amining  the  linear  and  quadratic  dependence  of  the  noise 
floor  with  respect  to  power.  The  linear  component  is  due 
to  shot-noise,  while  the  quadratic  variance  is  due  to  the 
intensity  noise  fluctuations  (see  equation  (A37)).  The  re¬ 
sults  are  shown  in  figure  7a  and  c.  The  noise  floor  was 
seen  to  only  increase  linearly  with  laser  power,  confirm¬ 
ing  the  absence  of  intensity  noise  at  the  frequencies  of 
interest. 

A  second  measurement  is  done  to  characterize  the 
phase  noise  properties  of  the  system.  By  sending  the 
laser  through  a  Mach-Zehnder  interferometer  (MZI)  with 
transfer  function  I(t)  =  /o(l  +  sin(27rcd/u;FSR)),  the  in¬ 
tensity  of  the  transmitted  light  will  contain  fluctuations 
related  to  the  frequency  fluctuations  of  the  light  (see  fig¬ 
ure  7).  The  free  spectral  range  (FSR)  of  the  MZI  is 
cdFSR/27T  =  115  MHz.  For  a  real  detector,  and  assum¬ 
ing  uj  <  g^fsr?  we  arrive  at 

Sll(cu)  =  S'dark(^)  +  |0TC>|2  +  ?7de  J  2 ~ ~  )  . 

Some  phase  noise  was  detected,  as  shown  in  figure  7b 
and  d  and  the  quadratic  dependence  of  the  PSD  on  sig¬ 
nal  power.  The  spectral  densities  show  a  roll-off  due  to 
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Figure  A7:  Laser  noise  characterization,  a,  We  measure  the  power  spectral  density  (PSD)  of  our  laser  for  several  powers, 
normalize  to  and  then  subtract  it  from  the  dark  current  of  the  detector.  The  same  measurement  is  performed  using  a  Mach- 
Zehnder  interferometer  locked  at  half  of  the  fringe  amplitude  in  order  to  convert  any  frequency  noise  to  intensity  noise  to  allow 
detection  and  is  shown  in  b.  c,  Plot  of  the  mean  value  of  the  PSD  around  the  mechanical  frequency  cum  from  the  measurement 
done  in  a  as  a  function  of  power.  The  good  linear  fit  (red  line)  indicates  that  no  intensity  noise  is  present,  d,  Mean  PSD  of 
the  measurement  in  b.  The  quadratic  fit  (red  line)  shows  that  phase  noise  is  indeed  present  (see  text  for  more  details). 


the  FSR  of  the  MZI.  It  was  found  that  in  the  frequency 
range  of  interest,  1  MHz  <  uj/2tt  <  40  MHz,  the  fre¬ 
quency  noise  spectral  density,  =  uj2 S </></> (co>) ,  is  flat, 

and  roughly  equal  to  3  —  6  x  103  rad2  -  Hz,  in  agreement 
with  previous  characterization  of  the  same  lasers  at  higher 
frequencies  [41]. 


C.  Linearity  of  detector  with  local  oscillator  power 


Having  characterized  the  laser  with  an  independent 
setup,  we  try  to  understand  the  properties  of  the  mea¬ 
surement  system.  Our  first  measurement  is  designed  to 
characterize  the  linearity  of  the  detector  and  amplifier. 
With  /dc  =  0,  and  no  signal  in  the  signal  arm  of  the 
BHD,  we  expect  the  system  to  faithfully  reproduce  the 
relation  (35)  showing  a  linear  relationship  between  local 
oscillator  power  and  the  detected  signal  vacuum  fluctua¬ 
tion  (shot-noise)  noise  level.  It  is  observed  that  the  mean 
value  of  the  PSDs  linearly  depend  on  the  input  power  as 
expected  and  shown  in  figure  8.  This  indicates  that  our 
detector  (and  its  amplifier)  are  in  fact  linear.  The  red 
line  is  a  linear  fit,  with  a  coefficient  of  determination  of 
R2  =  0.9999.  Although  we  already  confirmed  that  no 
measurable  amount  of  intensity  noise  is  present  (cf.  fig¬ 
ure  7),  in  the  case  we  would  have  an  appreciable  amount 
of  noise  this  measurement  would  show  that  it  is  smaller 
than  the  CMRR. 


a  b 


Frequency  (MHz)  Local  Oscillator  Power  (mW) 

Figure  A8:  Noise  level  versus  power,  a,  Electronic  noise 
power  spectral  densities  from  the  balanced  homodyne  detector 
at  different  local  oscillator  powers  (under  a  balanced  condi¬ 
tion).  The  red  trace  corresponds  to  the  electronic  noise  floor 
with  zero  local  oscillator  power,  i.e.  the  dark  current,  b,  Mean 
value  of  the  power  spectral  densities  shown  in  a  as  a  function 
of  local  oscillator  power.  In  this  plot  the  electronic  noise  or 
dark  current  contribution  (0.12  pW/Hz,  shown  by  the  dashed 
black  line)  is  subtracted.  The  red  line  is  a  linear  fit,  which  has 
a  coefficient  of  determination  R2  =  0.9999.  The  local  oscilla¬ 
tor  power  used  in  the  experiment  presented  in  the  main  text 
corresponds  to  3.0  mW. 
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D.  Detected  noise  level  with  unbalancing 


in  a  modification  on  the  order  of  one  percent. 
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Figure  A9:  Amplifier  gain,  a,  Shown  are  the  power  spec¬ 
tral  densities  (PSDs)  of  the  local  oscillator  as  a  function  of 
the  balancing  of  the  optical  power  in  the  two  paths  of  the  ho¬ 
modyne  detector.  Each  trace  represents  a  different  ratio  of 
power  in  each  path.  These  traces  were  taken  with  a  local  os¬ 
cillator  power  of  3.0  mW,  as  used  in  the  experiment,  b,  The 
mean  value  of  the  PSDs  normalized  to  the  perfectly  balanced 
PSD  are  shown  as  a  function  of  the  difference  voltage  on  the 
two  photodiodes  in  the  homodyne  detector,  where  zero  volt¬ 
age  represents  perfect  balancing.  The  green  line  is  a  linear  fit 
to  the  data,  while  the  black  curve  is  a  quadratic  model,  which 
describes  any  classical  intensity  noise  that  could  cause  the  dif¬ 
ference  in  the  level  of  the  PSDs.  The  red  curve  is  the  sum 
of  the  two.  The  change  in  PSD  with  homodyne  unbalancing 
can  be  fully  explained  by  the  small  signal  gain  weakly  depen¬ 
dent  on  the  detector  unbalancing  (linear  fit)  and  no  classical 
intensity  noise  (as  previously  determined). 


A  second  measurement  with  vacuum  input  on  the  signal 
is  done  to  understand  how  the  amplifier  in  the  homodyne 
detector  depends  on  the  DC  level  of  the  electronic  sig¬ 
nal  after  the  photocurrent  is  subtracted.  Here  we  use  the 
variable  coupler  to  change  the  splitting  ratio  and  cause  an 
imbalance  between  the  optical  power  levels  in  the  arms. 
The  detected  noise  floors  are  shown  in  figure  9a,  and  the 
mean  detected  PSDs  are  shown  in  figure  9b,  normalized 
to  the  shot-noise  level.  We  find  that  at  larger  Wc  (re¬ 
lated  to  Jdc  linearly),  there  is  a  very  small  (<  2%)  drop 
in  the  gain  of  the  detector.  Using  a  linear  fit,  we  ex¬ 
tract  an  adjustment  to  the  gain  vs.  output  DC  current. 
This  means  that  for  a  measured  noise  power  spectral  den¬ 
sity  Smeas(uj)  taken  at  a  DC  voltage  Vdc,  we  estimate 
that  the  actual  PSD,  compensating  for  modified  gain,  is 
S(u)  =  (1  +  Udc/(— 0.0096))_15'meas(^)-  This  modifica¬ 
tion  is  used  from  here  on,  and  only  reduces  the  amount  of 
squeezing  we  observe,  as  the  quadratures  with  squeezing 
are  always  at  positive  voltages.  Additionally,  the  largest 
DC  voltages  we  work  at  are  roughly  ±1  V,  which  results 


Figure  A10:  Detuned  noise.  The  laser  is  detuned  with  re¬ 
spect  to  the  cavity  resonance  by  1  nm  and  spectra  are  taken 
using  the  homo  dyne  detector  over  a  range  of  phase  angles, 
with  a  local  oscillator  power  of  3.0  mW.  This  lets  us  estimate 
the  amount  of  additional  intensity  noise  we  might  acquire  in 
our  optical  signal  train,  a,  The  dotted  blue  line  shows  the 
amplitude  of  the  interference  of  the  signal  and  local  oscilla¬ 
tor  as  a  function  of  time.  We  lock  at  several  relative  phases 
(color-coded  from  green  to  red  in  a  and  b)  and  plot  the  asso¬ 
ciated  normalized  power  spectral  densities  (PSD)  relative  to 
shot-noise  in  b.  For  every  second  measurement  we  switch  the 
signal  beam  off  to  obtain  the  shot-noise  level  (blue  traces). 
The  maximum  difference  in  the  noise  level  is  around  0.5%. 


E.  Estimating  added  noise  in  the  optical  train 


In  our  third  measurement,  we  reflect  the  laser  light  off 
the  end  mirror  of  the  waveguide  coupler  (detuned  by  1  nm 
from  the  cavity),  and  measure  the  detected  noise  level  as 
a  function  of  #,  the  phase  difference  between  local  oscil¬ 
lator  and  signal.  This  measurement  is  sensitive  to  both 
the  conversion  of  phase  noise  to  intensity  noise  through 
dispersion  in  the  optical  train,  and  added  noise  due  to  ad¬ 
ditional  noise  processes  in  the  optical  train  such  as  guided 
acoustic- wave  Brillouin  scattering  (GAWBS)  [46],  which 
could  cause  uncorrelated  noise  in  the  local  oscillator  and 
signal  arms  (see  Eq.  (A35)).  The  results  of  this  measure¬ 
ment  are  shown  in  figure  10a  and  b.  The  first  figure  shows 
the  DC  interference  signal  between  the  local  oscillator  and 
signal  used  in  the  measurement.  The  LO  power  used  for 
this  experiment  was  the  same  as  for  the  actual  squeezing 
data,  and  the  signal  level  used  is  on  the  same  order  as 
used  for  the  highest  power  measurements,  as  is  evident 
from  the  swing  of  about  1.4  V  in  the  DC  interference  sig¬ 
nal.  The  highest  DC  swing  observed  in  the  experiment 
was  1.6  V.  The  second  figure  shows  the  normalized  (to 
shot-noise)  power  spectral  density  where  an  increase  of  at 
most  0.5%  is  observed,  indicating  these  sources  of  noise 
do  not  contribute  in  our  experiment. 
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F.  The  effect  of  laser  phase  noise 

Using  the  measured  value  for  the  spectral  density  of 
phase  and  frequency  fluctuations  from  section  IV  B,  the 
effect  of  laser  technical  noise  on  the  detected  squeezing 
spectra  can  be  calculated.  Following  the  derivation  in 
section  I C  and  taking  the  classical  noise  component  of  the 
field  input  to  the  cavity  to  be  =  ia-iri(j)(uo)  (with 

a  corresponding  LO  phase  noise  of  a^o  (co)  = 

we  arrive  at  an  expression  for  the  output  noise  due  to 

input  phase  noise  from  the  cavity: 

«iut(w)  =  *«in(l  +  M  (w)  -  A2(lo))4>(uj).  (A38) 

Without  optomechanical  interaction  ( G  =  0)  we  find 
Ai(uj)  =  — ft/(i( A  —  uS)  +  ft/2),  and  A2{uS)  =  0.  We  cal¬ 
culate  the  expression  for  the  current  noise  due  to  laser 
phase  noise  using  this  expression: 

/(N)(w)  =  alo  aiutM +«LO[aiut)(-w)]* 

+aoutaLO  (W)  +  aout[«LO  (_W)]* 

=  F(u)4>(u)  (A39) 

with  F(uj)  =  i\ctLO&in\ [e~l0(r(uj)  —  r(0))  +  eie(r(— uS)  — 
r(0))*].  The  PSD  of  the  photocurrent  due  to  phase  noise 
is  found  to  be 

S^\w)  =  \F(w)\^{w).  (A40) 

For  a  system  with  no  dispersion,  r(uj)  =  const.,  it  can  be 
easily  shown  that  F(cu)  =  0  as  expected.  For  an  over¬ 
coupled  cavity  with  no  optomechanical  coupling,  r(ui)  = 
1  —  ft/(i( A  —  uS)  +  ft/2),  so  r{uS)  —  r( 0)  ~  4 za;/ft,  and  we 
have  F(cj)  =  8i|<aLO<ain|  sin(0)(u;/ft).  The  uo  dependence 
of  F(<jS)  means  that  a  flat  frequency  fluctuation  spectrum 
(S^  oc  u;-2,  as  we  observe)  adds  a  flat  noise  floor  to  the 
detected  signal. 

Finally  we  note  that  phase  noise  on  the  laser  can  drive 
the  mechanical  motion  and  cause  heating.  This  effect  is 
negligible  since  we  are  tuned  near  resonance,  where  only 
the  intensity  fluctuations  affect  the  mechanics,  and  our 
cavity  has  a  very  large  linewidth  ft. 

G.  Phenomenological  dispersive  noise  model:  the 
effect  of  structural  damping 

Mechanical  damping  of  resonators  and  the  associated 
fluctuations  from  coupling  to  the  thermal  bath  has  long 
been  considered  as  an  impediment  to  measuring  weak 
forces  in  gravitational  wave  detectors  [33,  34,  47-49].  In 
these  studies  the  effect  of  the  bath  has  often  been  en¬ 
capsulated  in  a  parameter  4/(u;),  representing  the  lag 
angle  in  the  response  of  the  material  to  a  force.  This 
lag  angle  is  the  complex  part  of  the  spring  constant: 
F  =  —k(l-\-i^f(u))x.  The  quality  factor  of  the  res¬ 
onator  is  given  by  the  narrow-band  properties  of  the  lag 
angle  and  its  value  at  the  mechanical  resonance  frequency, 


Q  =  4>(c Jm0)_1-  We  are  interested  in  the  wideband  prop¬ 
erties  of  4>(cj),  since  the  spectral  properties  of  the  thermal 
fluctuations  are  related  to  the  spectrum  4/(c<;),  following 
the  fluctuation-dissipation  theorem. 

In  the  case  of  our  experiments,  we  observed  noise  floors 
for  Sxx  following  a  cj-1  power  law  on  the  low  frequency 
end.  This  sort  of  noise  power  law  corresponds  to  a  flat 
spectrum  for  the  lag  angle  4/(tj)  =  const,  over  the  fre¬ 
quency  range  of  interest.  Unlike  viscous  damping  which 
can  be  simply  shown  to  have  4/(cj)  oc  c o  (since  the  force 
is  proportional  to  velocity),  a  lag  angle  constant  in  fre¬ 
quency  lacks  a  simple  physical  explanation,  though  it  is 
ubiquitous  in  many  types  of  mechanical  resonators  and 
commonly  called  “structural  damping”  [34]. 

In  the  input-output  formalism  outlined  in  section  I  we 
model  this  type  of  noise  by  taking  the  mechanical  damp¬ 
ing  rate  71  to  be  spectrally  flat,  and  using  frequency 
dependent  bath  correlation  functions  (b[n(uj  )&!>')>  = 
(n(u>)  +  l)S(w  +  a/),  ( bln(Lo)bin(u> '))  =  n(u>)S(uj  +  w'), 
{bjn((jj)bln(uj'))  =  0,  and  (&in(w)&in(a/))  =  0.  This  con- 
stitutes  our  single-mode  thermal  noise  model. 

In  any  real  optomechanical  system,  a  family  of  mechan¬ 
ical  modes  couples  to  the  optical  resonance.  In  the  modal 
picture  which  we  use  here,  each  of  these  mechanical  reso¬ 
nances  can  be  thought  to  add  to  the  detected  noise  floor 
with  its  contribution  scaling  at  the  low-frequency  end  as 
l J-1.  The  contribution  of  each  mode  is  proportional  to 
the  bath  temperature,  #Qk,  7^,  and  cj“2k.  We  lump 
all  of  these  contributions  into  a  single  effective  mechani¬ 
cal  resonance,  with  its  properties  (not  all  independent) 
determined  by  fitting  to  the  low  frequency  end  of  the 
noise  floor.  This  mechanical  resonator  is  modeled  with 
a  mechanical  frequency  ujm/ 2tt  =  50  MHz  (so  we  oper¬ 
ate  in  the  low  frequency  tail) ,  a  mechanical  quality  factor 
Qm  =  100,  and  a  total  coupling  rate  of  go/2ir  =  100  kHz. 
We  found  that  this  model  reproduced  the  magnitude  and 
phase  (the  quadrature  in  which  the  noise  is  detected)  of 
the  l j-1  noise  well,  if  an  additional  intracavity  photon- 
dependent  heating  of  Co  =  3.2  x  10-4  K/photons  is  as¬ 
sumed.  These  background  noise  floors  are  plotted  in  fig¬ 
ure  11.  This  cavity  heating  rate  leads  to  the  effective  bath 
temperature  to  nearly  double  at  the  highest  input  powers, 
going  from  16  K  to  over  30  K.  This  amount  of  heating  is 
in  line  with  what  we  expect  from  thin-film  photonic  crys¬ 
tals  we  have  fabricated  in  the  past  operating  in  the  same 
cryostat  [24]. 

H.  Phenomonological  absorptive  noise  model 

In  addition  to  the  noise  in  the  quadrature  of  the  me¬ 
chanical  motion  (which  arises  from  fluctuations  in  the 
cavity  frequency  cj0,  and  we  suspect  is  mechanical  in  ori¬ 
gin),  we  observed  a  significant  amount  of  noise  in  the 
opposite  quadrature,  which  can  be  interpreted  to  arise 
from  fluctuations  of  the  cavity  decay  rate  ft.  Additionally, 
we  observed  a  different  noise  floor  power  law  (cj-1/2)  for 
this  noise,  which  may  rule  out  an  optomechanical  origin. 
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Figure  All:  Noise  model  and  experimental  results.  The 

complete  noise  model,  and  constituent  components,  are  plot¬ 
ted  and  compared  to  the  experimental,  shot-noise-subtracted 
power  spectral  density  (PSD)  for  a  quadrature  sensitive  to  the 
mechanical  motion  (solid  curves)  and  an  insensitive  quadra¬ 
ture  (dashed  curves).  The  black  lines  are  the  experimental 
PSDs.  The  red  lines  represent  the  full  noise  model  includ¬ 
ing  contributions  from  a  single  mechanical  mode  (blue  line), 
phase  noise  of  the  laser  (brown  line),  and  the  extra  thermal 
noise  (green  line)  as  described  in  section  IV  G.  The  deviation 
between  the  modeled  and  experimental  data  predominantly 
results  from  additional  mechanical  modes. 


Figure  A12:  Power  spectral  densities  (PSD)  of  noise 
contributions  with  varying  powers.  The  complete  noise 
model  along  with  its  constituent  components  and  experimen¬ 
tal  data  are  shown  for  varying  optical  powers.  All  curves  have 
been  normalized  to  the  shot-noise  level.  The  experimental 
data  are  shown  in  black  with  the  full  noise  model  in  red  con¬ 
sisting  of  the  single  mechanical  mode  (dashed  blue),  phase 
noise  (dashed  brown) ,  and  extra  thermal  noise  (dashed  green) . 
The  optical  power  scaling  is  represented  by  the  transparency 
of  the  individual  curves  with  curves  becoming  less  transparent 
with  increasing  optical  power. 


The  power  law  scaling  agreed  with  thermorefractive  noise 
studied  extensively  in  the  context  of  gravitational  wave 
detection  [48],  microspheres  [50],  and  microtoroids  [51], 
but  it  is  expected  that  thermorefractive  coupling  is  pre¬ 
dominantly  in  the  same  quadrature  as  the  mechanical 
noise,  which  is  not  observed  here.  Also,  if  thermorefrac¬ 
tive,  the  noise  should  show  strong  variation  with  temper¬ 
ature  through  both  a  quadratic  temperature  scaling  (T2) 
and  an  extremely  steep  variation  of  dn / dT  in  the  temper¬ 
ature  range  of  16  K  to  30  K  [52],  which  was  not  observed. 
At  this  point,  we  have  no  noise  model  to  explain  the  ob¬ 
served  fluctuations,  and  the  origin  of  this  noise  will  be  the 
subject  of  further  investigation  to  be  presented  at  a  later 
time.  A  phenomenological  noise  model  was  instead  used, 
where  fluctuations  in  the  cavity  linewidth  proportional 
to  the  intracavity  power  with  aw-1/2  noise  spectrum  are 
assumed. 


I.  Comparing  measured  spectra  to  theoretically 
predicted  spectra 

Our  spectrum  analyzer  (Tektronix  RSA3408B)  oper¬ 
ates  by  taking  Fourier  transforms  of  a  time  domain  signal. 
By  windowing  a  short  time  sample,  and  calculating  its 


energy  spectrum,  a  power  spectral  density  is  constructed. 
The  size  of  the  window  in  the  time-domain  affects  the 
resolution  bandwidth,  and  is  well  known  in  signal  pro¬ 
cessing,  multiplication  by  a  Gaussian  window  of  length  r 
is  equivalent  to  convolution  of  the  frequency  domain  sig¬ 
nal  by  a  Gaussian  with  width  proportional  to  r_1.  All  of 
the  measured  data,  except  that  presented  in  section  III  C 
was  taken  with  a  40  MHz  window  and  300  kHz  resolution 
bandwidth.  The  spectra  contain  501  points  spaced  by 
80  kHz  in  the  frequency  domain.  Additionally,  for  a  few 
data  sets  we  took  narrow  band  spectra  (down  to  100  Hz 
resolution  bandwidth)  and  found  that  the  results  agreed 
over  the  regions  where  squeezing  was  observed.  The  the¬ 
ory  was  calculated  at  100  times  finer  resolution  than  the 
sampled  data  (with  50,000  points),  and  was  then  down- 
sampled  after  a  Gaussian  convolution  step  simulating  the 
operation  of  the  spectrum  analyzer.  This  only  effects  the 
size  of  the  mechanical  peak,  and  has  no  effect  on  the 
frequency  ranges  where  we  see  sub-shot-noise  fluctuation 
spectra.  For  the  thermometry  data  in  section  IIIC,  since 
we  are  interested  in  the  mechanical  linewidths  and  areas, 
the  span  was  always  chosen  to  be  the  minimum  allowable 
by  the  RSA,  which  is  twice  as  large  as  the  linewidth. 
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Figure  A13:  Detected  noise  power  at  a  given  frequency  vs.  the  lock  angle.  In  these  plots,  a  series  of  traces  is  shown 
of  the  detected  noise  level  at  a  given  frequency  (with  resolution  bandwidth  of  300  kHz)  as  a  function  of  the  locked  phase  6>iock- 
The  grey  points  are  the  measured  data  points.  The  solid  lines  are  the  results  of  models  detailed  in  this  text,  and  the  dashed 
lines  represent  different  components  of  noise  present  in  each  model.  The  red  line  shows  the  full  noise  model,  containing  the 
transduced  thermal  brownian  motion  from  the  studied  mode,  the  noise  due  to  structural  damping  present  in  the  system,  the 
phase  noise,  and  the  phenomological  out-of-quadrature  noise.  The  green  line  is  for  a  model  considering  all  the  same  noise 
contributions,  except  the  phenomological  component.  A  model  considering  a  system  without  any  thermal  noise  is  shown  in 
blue.  With  no  thermal  force  on  the  mechanical  systems,  the  detected  signal  in  this  case  can  be  attributed  to  radiation  pressure 
shot-noise  heating.  The  shot-noise  level  is  denoted  by  a  grey  line.  The  contribution  due  to  thermal  motion  of  the  mode  of 
interest  is  shown  by  the  dashed  blue  line.  The  noise  contributions  due  to  phase  noise  and  structural  damping  are  much  smaller 
and  shown  by  the  green  and  red  dashed  lines,  respectively. 


V.  SUMMARY  OF  NOISE  MODEL  spectra  shown  in  the  main  text. 

In  Table  I  we  present  a  summary  of  the  parameters 
used  in  the  theoretical  model  for  the  wideband  squeezing 
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Figure  A14:  Close-up  of  detected  noise  power  at  a  given  frequency  vs.  the  lock  angle.  This  close-up  shows  regions 
of  squeezing,  and  the  colors  are  the  same  as  in  Figure  13. 


VI.  MATHEMATICAL  DEFINITIONS 


We  present  here  the  notational  conventions  used 
throughout  this  work  for  reference.  The  Fourier  and  in¬ 
verse  Fourier  transforms  of  operator  A(t)  are  defined  as 


1  r°°  .  . 

A(t)  =  -=  /  duo  e~lujtA(u)  and 

i  r°° 

4M  =  ^=  dt  elutA(t),  (A41) 

V  47T  J  —  (X) 


respectively.  The  Hermitian  conjugate  of  operator  A(t) 
is  given  by  A\t)  which  has  the  Fourier  transform  A^(u). 
This  is  related  to  A(uS) 

(^(w))1  =it(-w).  (A42) 

In  the  derivations  presented  here,  we  typically  express  a 
given  operator  in  terms  of  the  “input”  bath  operators. 
Expected  values  are  then  defined  as  (. A(t ))  =  Tr[pinA(t)], 
where  p-in  is  the  mixed  state  describing  the  bath.  Spectral 
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TABLE  I:  Model  Parameters 


Symbol 

Name 

Value 

Measurement 

^optical 

Optical  quality  factor 

5.7  x  104 

Low-power  optical  spec¬ 
troscopy  with  waveme- 

ter.  Ill  B 

V* 

Cavity- waveguide  coupling  efficiency 

0.55 

Low-power  optical  spec¬ 
troscopy  with  wavemeter. 
Verified  phase  response  with 
VNA  to  distinguish  from 
under-coupling.  IIIB 

7i/2  7T 

Mechanical  linewidth 

172  Hz 

Linewidth  measurement  vs. 
laser  detuning  in  thermom¬ 
etry  measurement  (see  sec¬ 
tion  III  C). 

go/2-K 

Optomechanical  coupling  rate 

750  kHz 

Linewidth  and  mechanical 
frequency  measurement  vs. 
laser  detuning  in  thermom¬ 
etry  measurement  (see  sec¬ 
tion  III  C). 

rriO 

1b 

Bath  temperature 

16  K 

Calibrated  areas  in  ther¬ 
mometry  measurement  (see 
section  III  C). 

CO 

Heating  by  optical  absorption 

3.2  x  10-4  K/photon 

Rise  of  uj-1  noise  floor 
with  optical  power  (see  sec¬ 
tion  IV  G) .  The  cavity  tem¬ 
perature  according  to  this 
model  rises  from  16  K  to 
roughly  30  K  at  the  highest 
powers. 

Sujuj 

Frequency  noise  spectral  density 

6  x  103  rad2  Hz 

Frequency  noise  measure¬ 
ment  with  Mach-Zehnder 
Interferometer  (see  sec¬ 
tion  IV  B). 

A 

Laser  detuning  (red  laser  is  positive) 

(0.044  ±  0.006)k 

The  intensity  of  the  reflected 
light  is  used  to  initially  set 
the  detuning.  For  a  more 
accurate  determination,  the 
value  of  A  minimizing  the 
detected  signal  for  the  ob¬ 
served  0*ock  is  found  (see  Sec¬ 
tion  II C). 

$lock 

lock  angle 

varies 

The  lock  point  (as  in  fig¬ 
ure  3)  is  used  to  find  the 
phase  angle  between  the 
light  reflected  from  the  cav¬ 
ity  and  the  local  oscillator. 

^lock 

critical  lock  angle 

varies 

This  is  the  lock  angle  were 
no  mechanical  signal  is  de¬ 
tected.  It  is  found  by  looking 
at  the  area  of  the  mechanical 
mode  as  a  function  of  0\oc k 
(see  Section  II C). 

densities  are  found  by  taking  the  Fourier  transform  of  the 
auto-correlation 


For  a  real  operator,  Saa{u)  =  Saa(~u).  For  non-real 
operators,  the  spectral  density  can  be  symmetrized  to 


/oo 

dr  el“T{A]{t  +  r)A(t)). 

-OO 


(A43) 
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Saa(w): 


We  use  the  convention 


Saa(A  =  \(SAa{A  +  Saa{~A)-  (A44) 

These  spectral  densities  can  also  be  calculated  from  the 
frequency  domain  operators: 

/oo 

dJ  (it(w)i(a/)).  (A45) 


X(ej)  =  aje-ie  +  a]  eie .  (A46) 

to  define  a  measured  quadrature  of  the  field.  This  defini¬ 
tion  corresponds  to  having  a  phase  difference  of  0  between 
input  light  and  local  oscillator. 
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